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CHAPTER I 


INTRODUCTION 

The ever-incrpasing demand for additional communications capacity 
has led system design engineers to increase operating frequencies 
higher into the millimeter wave frequency band. Also, to effectively 
double communications capacity for a given satellite communications 
system, frequency reuse systems employing orthogonal polarizations have 
been suggested [1]. However, for operating frequencies above 10 GHz 
rain attenuation and depolarization become more severe because of the 
in-reased signal scattering and absorption by raindrops. Before future 
coranuni cations systems can be designed for maximum reliability and 
economy, a clearer understanding of these weather induced phenomena is 
necessary. Presently an experimental data base is being collected to 
describe the effects of rain on terrestrial and satellite communications 
links. Measured signal attenuation, orthogonal channel isolation, and 
phase shift data are being correlated with weather data to aid the 
design engineer in the prediction of rain effects on miUimeter wave 
propagation. Since this collection of data is relatively expensive, 
especially for satellite communications links, theoretical models have 
been developed to predict weather effects on communications link 
performance. This report investigates the ability of a theoretical 
model to account for the inhomogeneous vain aonditicns present on a 
satellite communications link and thus accurately predict the effects 
of weather on system performance. 

Theoretical modeling of the scattering properties of an ensemble 
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of particles was considered by Gustav Mie as early as 1908 [2]. The 
present development of theoretical modeling as related to a rain-filled 
space and the restrictions associated with the ensemble of raindrops 
will be discussed in Chapter II. In Chapter III a new deterministic 
theoretical model is presented with an associated rain propagation 
prediction program that can model an inhomogeneous rain and accurately 
predict the effects of precipitation on a satellite communications link 
for a variety of frequencies, elevation angles and locations. Finally 
in Chapter IV, the most recent multiple frequency data from various 
depolarization experiments are presented as verification of the new 
model . 

The new model presented in Chapter III, although more general than 
previous models, is theoretically equivalent to each of the classical 
models presented in Chapter II under the assumptions of the respective 
model. By taking a slightly different approach in the theoretical 
derivation of the new model, the scattering properties of an arbitrary 
inhomogeneous ensemble of raindrops and ice crystals can be determined 
and the validity of other theoretical models can be verified. Unique 
aspects of the new model are: 1) attenuation and depolarization are 

computed directly in terms of the scattered fields from an ensemble of 
particles containing an arbitrary mixture of particle type (rain or 
ice), particle shape, particle size, particle density and particle 
orientation, 2) an inhomogeneous mixture of particles and the presence 
of a melting layer along the propagation path is accounted for directly 
in determining the scattering fields exiting the rain medium, 3) a 
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frequency independent synthetic storm algorithm models the nonuniform 
rain rates present on a satellite communications link, 4) antenna-wave 
interaction is considered, and 5) signal attenuation, channel isolation 
and phase as a function of ground rain rate can be determined for a 
variety of site locations from a knowledge of the operation frequency, 
system antenna polarization parameters, and the system elevation angle. 



CHAPTER II 


REVIEW OF CURRENT THEORETICAL PROPAGATION MODELS 

Although Gustav Mie in 1908 [2] was responsible for the prelimi- 
nary foundation for theoretical work relating to the attenuation and 
depolarization of electromagnetic waves propagating through rain, it 
was not until 1960 that Oguchi [3] first considered the depolarization 
of an incident field by an oblate spheroidal raindrop. Oguchi obtained 
the first-order change in the scattered field from a single oblate 
raindrop by expanding the scattered field directly in a series of 
spherical vector wave functions involving powers of the raindrop 
eccentricity. The single oblate raindrop scattering coefficients 
provided the foundation for all theoretical rain depolarization models 
that followed. The scattering coefficients of Oguchi were later veri- 
fied using different theoretical methods by Morrison and Cross of Bell 
Labs in 1974 [4] and Uzunoglu, Evans, and Holt of the University of 
Essex, England in 1977 [5]. 

Using the theory of van de Hulst [6], Oguchi extended the single 
drop solution to an ensemble of identical, equioriented oblate sphe- 
roidal raindrops. By decomposing the incident field along the ortho- 
gonal axes defined by the principal planes of the canted oblate rain- 
drop, Oguchi determined an effective complex refractive index of the 
rain medium. The application of the work of Oguchi to the problem of 
depolarization as related to a communications system was first made by 
Thomas of Bell Labs in 1971 [7]. He showed that a difference In 
attenuation for waves polarized along the major and minor axes of an 
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oblate raindrop would lead to depolarization of any wave not polarized 
parallel to one of the principal axes. However, Thomas neglected the 
effects of differential phase on the depolarization of an electric 
field. The assumptions of the theoretical model of Thomas are given 
below: 

1) Uniform rainfall over the entire path, 

2) Laws and Parsons [8] raindrop-size distribution, 

3) Single scattering {which is to say that each raindrop 
is illuminated by an incident plane wave and wave 
interaction between drops is negligible), 

4) 100% oblate spheroidal equioriented raindrops at an 
effective canting angle B. 

The theoretical models of Watson (1973) [9] and Chu (1974) [10] also 
employ the assumptions of Thomas; however, these models consider the 
effects of differential attenuation and differential phase. The work 
of Watson and Chu, referred to as the classical propagation constant 
formulation, predicted rather well the rain effects experienced on a 
terrestrial communications link with nearly uniform rain conditions. 

Wiley [11] in 1973 checked the validity of the classical propaga- 
tion constant model. He derived the same results by directly summing 
the scattered fields located within the first Fresnel zone from an 
arbitrary equioriented ensemble of raindrops. (A drop size distribu- 
tion was not included.) Wiley also illustrated that theoretical models 
could represent an ensemble of raindrops that contained up to 60% 
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spherical raindrops and accurately predict measured depolarization data. 
This difference in opinion of the percent of oblate raindrops is 
complicated by the positions of meteorologists Jones [12] and Pruppacher 
and Fitter [13]. Using experimental vertical wind tunnel conditions 
and formulations relating surface tension to particle shape, Pruppacher 
and Fitter hypothesize that, except for only the largest drops, rain- 
drops are oblate spheroidal particles. Based on real rain measurements, 
Jones indicates that a real rain is composed of spherical, oblate 
spheroidal, prolate spheroidal and irregular shaped raindrops for all 
size classes greater than 0.9 mm spherical radius. Depending on the 
assumption of the percent of oblate raindrops present in an ensemble of 
raindrops, channel isolation can vary by a significant amount for a 
given value of attenuation as indicated in Fig. 2-1. (This figure was 
obtained from the results of Chapter III.) Since the work of Jones 
indicates the presence of shapes other than oblate and the effect of 
different raindrop shapes has a significant impact on the prediction of 
channel isolation, a theoretical propagation model should include a 
particle shape distribution. 

Although still assuming 100^ oblate spheroidal raindrops, Uzunoglu, 
Evans, and Holt [5] and Nowland, Olsen, and Shkarofsky [14] in 1977 
indicated that a distribution of canting angles should be included in a 
theoretical model. Since the classical propagation constant formula- 
tion is restricted to equioriented raindrops, a stochastic averaging 
process was implemented. 3y weighting the elements of the scattering 
matrix as seen in Eq. (2-1) with a Gaussian probability distribution 



ISDLflTlUi'i (DB) 




- 8 - 


function, a canting angles distribution can be modeled: 


2wDii(Q*) swD-j2(®.j) 


2 w D2 i(0^-) 
i 


I w DjjCe,) 




where 


D-J1 = dj^ cos^ ® ® 

Doo = dr sin^ 0 + d cos^ 0 (2-1) 

I'd n V 

°12 *^21 ^ ® ® 

and the symbols w, d^ 0, N, and L are defined as 

w = normalized Gaussian weighting function 

-jXL I N(a) f^^ r(a) 

d . = e ' * 

v,h 

6 = drop canting angle 

fy h ® principal plane single drop scattering coefficients 
N(?) = number of raindrops per unit volume in the a" to ? + 
da equi volumetric drop radius interval. 

L = rain extent length 


Because a stochastic averaging process was used, the above formulation 
ignores multiple scattering between arbitrarily oriented raindrops. 
However, the assumption was made that multiple scattering is negligi- 
ble. This assumption can be verified by the more general model 
presented in the following chapter because it accounts for the arbi- 
trarily oriented particles directly. 
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For satellite communications links the assumption of uniform rain- 
fall ever the entire rain extent is not a.aquate in describing the 
effects of precipitation on system performance. As peak rain rate 
increases, rain cells of higher rain rates decrease in size. Based on 
experimental attenuation measurements, researchers [15], [16], [17] 
have implemented an effective path length (L^) that models the nonuni- 
for rain rates present on a satellite conmuni cations link. However, 
it will be shown in Sec. 3.6 that this method is frequency dependent 
and overly restrictive. Therefore, a frequency independent synthetic 
storm algorithm will be presented. 

The formulations of Uzunoglu, Evans and Holt [5] and Howland, 
Olsen, and Shkarofsky [14] represent the current status in theoretical 
propagation models. However, these formulations do not consider the 
effects of an ensemble of arbitrary particles and an inhomogeneous rain 
medium. The theoretical model presented in the following chapter 
accounts for the scattering properties of an arbitrary inhomogeneous 
ensemble of particles without a significant increase in model complex- 
ity. 





CHAPTER III 


DEVELOPMENT OF THE THEORETICAL MODEL 

A theoretical rain propagation model should represent accurately 
the physical na,ture of a rain medium and allow a wide range of rain 
parameters; thus, the effect of precipitation on communications system 
performance can be determined accurately. Parameters such as canting 
angle, particle size, particle shape and particle type should be 
described with deterministic distributions based on existing physical 
knowledge of their behavior. The model also should be flexible so 
that these parameters can be adjusted to best describe the true 
physical nature of the rain medium. In this chapter a model is 
presented that provides this generalized format and henceforth will be 
referred to as the scattering modet. 

3.1 Electromagnetic Scattering by an Oblate Spheroidal Particle 

To develop a solid foundation for the derivation of the scattering 
model, the scattering properties of a single oblate raindrop will be 
considered first. The oblate raindr’op is of extreme importance in any 
model that predicts the attenuation, depolarization, and phase shift 
of an electric field. The oblate particle through its differential 
attenuation and differential phase properties provides the key mecha- 
nism in the depolarization process. Oguchi and Hosoya [18], Morrison 
and Cross [4], and Uzunoglu, Evans, and Holt [5] have published single 
drop scattering coefficients as a function of drop size for the princi- 
pal planes of the oblate drop. The drop size is represented by an 
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equi volumetric drop radius, and as the radius increases the drop 
becomes more oblate. The available scattering coefficients agree 
rather well and provide the basis for any theoretical model. 

An oblate spheroidal raindrop canted at an angle 9 with respect to 
an arbitrary x-y coordinate system is illustrated in Fig. 3.1-1. If 
the x-y coordinate system is aligned with true horizontal and vertical, 
6 is the drop canting angle. Following the notation of Oguchi, the 
h-v axis in Fig. 3.1-1 represents the orientation of the oblate rain- 
drop where the h and v axis are aligned with the drop major and minor 
axes, respectively. Another important parameter needed in the calcula- 
tion of the single drop scattering coefficients is the elevation angle. 
The elevation angle s is the angle between a plane parallel to the 
local horizontal and the direction of propagation represented by the 
vector K^. For a terrestrial consnuni cations link, this angle would be 
zero. A satellite link would have an elevation angle ranging from 0 to 
90 degrees. Knowing the equi volumetric drop radius as defined by 
Uzunoglu, Evans and Holt and given in Fig. 3.1-1, the elevation angle, 
and the frequency of propagation, the vertical and horizontal scatter- 
ing coefficients f^ and f^ can be determined. 

Any incident field can always be decomposed into orthogonal linear 
components along the principal axes of the oblate spheroidal raindrop. 
Thus, single drop scattering coefficients are given only for the 
principal axes of the drop. Single drop scattering coefficients are 
defined in several ways in the literature; the definition of Oguchi 
will be used as follows: 
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Figure 3.1-1. Arbitrary oblate spheroidal 

raindrop canted in an arbitrary 
x-y coordinate system. 
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f : ratio of the forward scattered electric field 

V 

component along the drop minor axis to the 
incident electric field component along the 
drop minor axis 

f^: ratio of the forward scattered electric field 
component along the drop major axis to the 
incident electric field component along the 
drop major axis . 

The incident electric field component along the drop major axis (h) 
experiences more attenuation and phase shift than the component along 
the minor axis (v). This effect causes the depolarization of incident 
fields that are not linearly polarized along one of the principal axes. 

The hv decomposition is overly restrictive. Wiley [11] has shown 
that a generalized scattering coefficient can be defined as a function 
of 9 and the single drop scattering coefficients. The generalized 
scattering coefficients allow the incident field to have any arbitrary 
angle with respect to the drop. These coefficients are easily derived 
(see Appendix 6.1.1) and are given below in terms of the arbitrary x-y 
coordinate system: 

^xx * ^v ® ^h ® 

fxy = " f|,) sin 9 cos e 

f = f 
yx xy 

fyy * f^ cos2 9 + f sin2 6 . 


(3.1-1) 
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The symbol f will be used to represent the general scattering coeffi- 

r H 

dent. 

Now the forward scattering effects of an oblate spheroidal rain- 
drop can be examined. The incident field on the drop can be written in 
the general form 

E’ = Ejj J + Ey p (3.1-2) 

• • 

for all polarization states where E and E are phasor quantities. The 

X y 

forward scattered electric field at point P in Fig. 3.1-1 is expressed 
by the relation 




(3.1-3) 


where 


E (P) = 
pq^ ^ 


" i 

Ep(P) 


-jkor 


pq 


(3.1-4) 


= the scattered electric field component with 
polarization q as a result of an incident field 
with polarization p, evaluated at point P 


and, 


f«« = the ratio of the forward scattered electric field 

pq 

with polarization q to the incident electric field 
with polarization p evaluated for a particular drop 
orientation, drop size, and frequency, 

-jk^r 

The - — factor of Eq. (3.1-4) describes the spherical wave behavior 
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of the scattered field between the drop and point P. Although the 
oblate drop does have back and side scattering properties, these are 
of no consequence in the prediction of rain effects on a point to 
point millimeter wave conmuni cations link. However, back-scattering 
coefficients are available in the literature [18] and can be used in 
place of the forward scattering coefficients in Eqs. (3.1-1) and 
(3.1-4) to obtain the back-scattered electric field. 

The forward scattering effects of an oblate raindrop can be 
described by Eq. (3.1-3) in terms of the drop orientation, drop size, 
and frequency. The generalized scattering coefficients of Wiley [11] 
also can be extended to <my particle having two axes of symmetry. 

With this fact, the scattering effects of a slab of arbitrary particles 
can be investigated. 

3.2 Electromagnetic Scattering by a Single Homogeneous Slab of 

Arbitrary Particles 

To predict accurately the effects of precipitation on a given 
communications link, the propagation path may be subdivided into small 
Intervals of thickness which are commonly called rain slabs. To 
describe the scattering properties of a slab, and hence the scattering 
properties of the entire propagation path, a physical understanding of 
the rain parameters 1s r,ecessary, Jones [12] has indicated that rain 
consists of particles with a variety of shapes. An average rain 
consists of spherical, oblate spheroidal, prolate spheroidal, and 
irregular shaped raindrops. Furthermore, it is not uncommon to find 
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that a heavy concentration of ice crystals at higher altitudes affects 
the performance of a satellite link. The particle shape, particle 
orientation, and particle density influence electromagnetic waves 
propagating through a rain medium. All particle shapes except for 
spherical contribute to the depolarization of the incident field. For 
this reason the electromagnetic scattering from a single homogeneous 
slab of arbiirary particles is extremely important to the evaluation 
of rain effects on millimeter wave propagation. 

In 1957 van de Hulst [6] addressed the problem of light scattering 
by a slab of particles. Although the scattering here is at high radio 
frequencies, van de Hulst's optical work will serve as a foundation for 
the development of the scattering model. Consider a homogeneous slab 
of particles aa meters thick and infinite in extent in the transverse 
plane as shown in Fig. 3.2-1. Since single drop scattering coeffi- 
cients are not readily available for irregular and prolate spheroidal 
raindrops, the assumption will be made that all scattering effects of 
the slab are produced by spherical raindrops, oblate spheroidal rain- 
drops, and ice crystals. However, if coefficients do become available 
for the irregular and prolate spheroidal raindrops, their inclusion in 
the theoretical development of the model would further refine the 
physical representation of the rain medium. With the exception of the 
work of Wiley [11], most researchers have assumed 100% oblate raindrops 
which is not physically accurate. Also, any model with the 100% oblate 
assumption inherently predicts more depolarization for a given rain 
rate and attenuation because non-depolarizing spherical scatterers are 
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Particle 



Figure 3.2-1. Hcanogeneous rain slab of 
arbitrary particles. 
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omitted. 

The incident field on the slab of particles is assumed to be a 
uniform plane wave propagating in a direction perpendicular to the slab. 
Also, the magnitude of the incident field is the same for all particles 
within the slab. The incident field, regardless of polarization, can 
be represented by its x and y components as in the previous section. 

The rain slab is assumed to be large-scale homogeneous with a mixture 
of particles having varying shapes, sizes, and orientations within the 
volume of the slab. The distribution functions describing these para- 
meters will be discussed in greater detail in Sec. 3.6. The purpose of 
this section is to determine the total electric field exiting the propa- 
gation medium. 

The total electric field at point P is the sum of the incident 
field and the forward scattered field both evaluated at point P. 

F(P) = e'’(P) + ?(P) . (3.2-1) 

The total forward scattered field at point P is the sum of the scat- 
tered fields from all the particles within the slab. The total scat- 
tered electric field with polarization q resulting from an incident 
field with polarization p evaluated at point P is expressed by 

. .. __ ' ^m^ 

^ ^ ? * (3.2-2) 

P m PQm 

All Particles 

The factor f (i‘, T) is the generalized complex forward scattering 
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coefficient for the m particle as found in Eq. (3.1-1) and represents 
the change in amplitude, phase, and polarization that the incident field 
experiences as a result of the particle. Note that the scattering 
coefficient is a function of particle size (?) and particle orientation 


(?) and must remain within the summation. The factor e 






/r, 


m 


represents the spherical nature of the scattered waves and the effect 

th 

of the location of the m particle on the phase of the scattered field 
at point P. 

With a few assumptions Eq. (3.2-2) can be simplified. Consider 
the geometry for the location of the m^" particle as shown in Fig. 

3.2-1. The relative phase difference corresponding to the path 


difference is all that need be considered in determining the 

phase of the scattered waves arriving at point P. If the assumption is 

made that then r„ - z„ and: 

m ^m m m 


= p2 + 2^ 

m ^m m 


p2, = r^ - z2 = (r^ - z^)(r_ + z^) 
m m m m m' ' m m^ 


^’^m “ ^m^ " + z 1 ' 

^ m *m^ 


r,. , 


(3.E-3) 


The assumption » p^^^ is reasonable because the scattered fields at 
point P are Influenced coherently only by the particles located within 
the first few Fresnel zones. For a satellite link, the radius of the 
first Fresnel zone corresponding to a rain extent less than 15 km is 
less than 20 meters. Using Eq. (3.2-3) in Eq. (3.2-2) yields 
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P|/ 


Epq(P) » eJ(P) j f„„ (T. 0) ® 


2 z, 


in 


m 


All Particles 




(3.2-4) 


Although the evaluation of the sum in Eq. (3.2-4) is a relatively simple 
procedure, several new symbols are needed to abbreviate the somewhat 
long and cumbersome summations. These symbols will be defined as they 
are needed. 

Let Np represent the total number of particles within the slab, 
and let Pe^, Pe^, and Pej represent the effective fraction of particles 
that are spherical raindrops, oblate raindrops, and ice crystals, 
respectively. Since there are three distinct particle shapes, Eq. 
(3,2-4) can be rewritten as 




Spherical Drops 

Pei Np 

Ice Crystals 


tn 


where 


PeoNp 

z 

m-^l 


pOBL.r 
pq 


(a, ^g, 


m 


Oblate Drops 


m 


(3,2-5) 




0 


g » ^ 






The superscripts SPH, OBL, and ICE denote the scattering coefficient 
for the respective particle type. 

Within each class of particle shapes there exists a distribution 
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of particle size and for the oblate raindrop and ice crystal classes 

SPH OBL 

there is also a distribution of particle orientation. Let N!j , k’ 

tf*p 

and ^ represent the number of spherical, oblate, a^d ice particles 
within a rain slab that are in the i^^ size class and k orientation 
class, respectively. The scattered field can now be expressed as a 
series of sunmations: 



I ! 


fjOBL 

rOBL/- T \ + 

„ , pq ' 1 ’ ^m 
m=t 


^OBL 

pq ^ n’ 
m= 1 


N 


ICE 




1 si Trr 


N 


ICE 

n,1 


+ I f^^^fa , 6i) 9 

, pq ^“n’ V ^m 
m=1 


N 


ICE 

1 


n 


ICE 


z f^^^(a, , 6 ) g + ••• + s* f^^^(a , e ) g 


m 


‘‘ ' s «_/ a_ 

_-j pq i w 


m=l 


pq ' n’ m* 


(3.2-6) 
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Note that the sumnatioii over the spherical drops is Independent of 
canting angle. Although Eq. (3.2-6) Is somewhat cumbersome, it Illus- 
trates the physical significance of the particular class of particles 
involved in the summation process. The assumption has been made that 
the canting angle of a given particle is independent of the equivolu- 
metric particle radius. Brussaard [19] has shown that this is not 
strictly true. As the radius and the terminal velocity of the falling 
particle increases, wind gradients have more effect on the orientation 
of the particle. However, the assumption of independence is reasonably 
accurate and it also leads to a relatively simple model. 

Since the scattering coefficients in Eq. (3.2-6) are constant 

within their respective summations, they can be taken out of the summa- 

Nl,k 

tion signs leaving the summation E g^^ in every term. This suumation 
is easily evaluated (see Appendix 6.1.2) and the result is a constant 
given below: 


N. 


I.k 
z 


2 g- = n. . , 
111=1 ” 


(3.2-7) 


The symbol n.. is the number of particles per unit volume in the i 
size class and the k orientation class. Equation (3.2-6) now 
becomes ; 


.th 




Ej(P) { n f H . 


+ 7 (I ) ) 

n pq' n' ' 
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OBL iOBL/— ^ \ 4. 


OBL .0BL,~ r 1 
"n,m ^pq V 


' "1.1 ^pq (®T ®1' * •“ 


„ICE ftCE,- T- > 
n_ T__ \ 3_ s Q_ / 
n,OT pq n m' 


. (3.2-8) 


Recognizing that the individual bracketed terms are double sums over 
the discretized particle size and particle orientation distributions, 
Eq. (3.2-8) can be further simplified. Let Ny(¥., ?j^) be the total 
number of particles per unit volume with particle size and particle 
orientation Then Eq. (3.2-8) can be rewritten as follows: 

M_. M,. 


Epq(P) ' E’(P) N^(a,.. q^)Fp^{a^. 6^) (3.2-9) 


where , 


fpqf^i- " '’*0 C<"i- 


(3.2-10) 


and, the symbols M. and M|^ represent the number of discrete intervals 
in the respective distribution functions. 

Although the scattered field expression has been simplified to a 
great extent, it is convenient to define an eff active scattering coef- 
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ficlent for the entire slab, Sp^, such that 

E® (P) = S„ ^{P) . (3.2-11) 

pq pq p 

From Eq. (3,2-9), the effective scattering coefficient Sp^ is expressed 
by the relation 



jXA& 


z 

i=l 



Ny(a^. 6|t)Fpq(a^» 



(3.2-12) 


= the ratio of the forward scattered electric field 

with polarization q to the incident electric field 

/ 

with polarization p representing the scattering 
properties of a thin rain slab with arbitrary particle 
shapes, sizes, and orientations. 


Using Eq. (3.1-1) in Eq, (3,2-12), the effective scattering coeffi- 
cients relative to the arbitrary x-y axes can be formulated as 




^xy ” ^yx 


Syy = 


+ cos2 e^CPeg f®*- + Pej 


"l ”k 

-jxw E E SuKPegCf^ - fj,)™*- 

1*1 ft*i 

+ Pei(^'y - ff^)’^*'^)sin cos 

”i ”k 

1 I I 

+ sin2 e,j(Pe(, + Pej 


(3.2-13) 
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The cross polarization scattering coefficients S and S represent 

xy yx 

the scattering of energy from one polarization component into its 
orthogonal component. This effect is directly related to the differ- 
ential phase and differential attenuation (f^ - f^) properties of the 
oblate raindrop and the ice crystal. The spherical raindrop on the 
other hand does not contribute to the depolarization process and when 
PeQ and Pej are zero, =s 0, and the rain slab does not 

depolarize the incident field. 

Using Eqs. (3.2-1) and (3.2-11), the total field at point P now 
Gan be expressed in terms of the effective scattering coefficients as 
f ol 1 ows : 

E,(P) = E>) . E^(P) . E’(P) 

(3.2-14) 

Ej,(P) = eJ(P) + E^(P) + E^(P) Syy . 

In matrix form, Eg. (3.2-14) becomes 


“ — 


- 



E>) 

yp) 



^yx 


- 


^xy 

(1 + S^y) 


E)(P) 


The scattering properties of a single homogeneous thin rain slab 
of arbitrary particles now have been presented. Also an effective 
scattering coefficient has been formulated that includes the effects of 
arbitrary particles, a particle size distribution, and a particle 
orientation distribution. So now with an understanding of the scat- 
tering properties of a thin rain slab, the effect of a homogeneous 
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rain on a conmuni cat ions link can be investigated. 

3.3 Electromagnetic Scattering by a Hoinogeneous Rain 

The purpose of the scattering model is to predict the physical 
changes that an incident electromagnetic wave with a given polarization 
undergoes when all or part of the propagation path is filled with rain. 
Now that the scattering properties of this path are understood on a 
microscopic level, the entire rain extent can be modeled by extending 
the work in the previous two sections. In this section the assumption 
will be made that the rain is homogeneous in the direction of propa- 
gation. However, this restriction will be relaxed in the next section 
when a model for a piecewise homogeneous rain is developed. 

Consider a propagation path with a homogeneous rain of extent n 
meters in the direction of propagation. In the previous section the 
total electric field was determined at point P for a single thin rain 
slab. Since the rain volume is homogeneous in the direction of propa- 
gation, the rain extent a can be modeled by a series of Mg thin rain 
slabs. If the incremental length aa of the thin rain slab is assumed 
to be one meter (aa = 1 m), is a large number equal to the homo- 
geneous rain extent length a. 

As a result of the assumption of a homogeneous rain in the direc- 
tion of propagation, the Mg slabs have identical scattering properties 
as described by Eg. (3.2-15). The effect of Mg rain slabs on the 
incident field can be represented by the multiplication of the Mg 
matrices as follows; 
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\h 


■(1 + S^) 


«s 

“•iH 

b. 


^xy 

0 + 


y 


(3,3-1) 


The scattering matrix [S] ® in Eq. (3.3-1) represents the changes 
that the incident field experiences as it propagates through a homo- 
geneous rain medium of thickness A. However, raising a matrix to the 

t*K 

power where H can be a large number (from 1000 to 10000 depending 
on the length i) is numerically awkward. This problem can be solved 
by the use of the Cayley-Hamilton theorem [20], The Cayley-Hamilton 
theorem allows the 2 x 2 matrix of Eq. (3.3-1) to be expressed in terms 
of the identity matrix and two unique constants as follows: 



= + a-j[S3 


(3.3-2) 


where. 


M M 

Oq “ ^2 ” ^1 ^2 ^^^^2 ~ ^ 


H M 

~ ^^2 * ^1 " ^1^ 


(3.3-3) 


(S +S ) / (S +S ) 

M.2 = * V <1 Sjy 


(1 + S )(1 + S ) . 
•'xx^' yy' 


Thus, the electric field exiting the rain medium now can be represented 

by 
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{a )} 

^ 0 r 


a-i S 
1 yx 




GCi S 
1 xy 


<“o ^ * VU 


L^y. 


(3.3-4) 


In this section the scattering properties of a homogeneous rain 
have been modeled to predict the effect of this rain on an arbitrary 
incident field. Equation (3.3-1) provides a physical insight into the 
scattering properties of the medium while Eq. (3.3-4} provides an effi- 
cient method for computing the electric field exiting a homogeneous 
rain. 


3.4 Electromagnetic Scattering by a Piecewise Homogeneous Rain 

The development of the scattering model as presented in Sec. 3.2 
considered a homogeneous mixture of particles with varying sizes, 
shapes, and orientations in a thin slab of rain. In the previous 
section a homogeneous rain was modeled as a series of identical thin 
rain slabs. In order to permit more accurate modeling of a real rain, 
the homogeneous rain restriction is removed in this section. 

By discretizing the rain extent on a large scale basis as shown 
in Fig. 3.4-1, a piecewise homogeneous rain can be modeled with M 
discrete rain oelte. Each rain cell can have arbitrary length and 
arbitrary rain conditions. The rain cells also can have different 
distributions of particle shape, particle size, and particle orienta- 
tion. The scattering properties of each rain cell can be represented 
by a scattering matrix [S.] where the effect of the i rain cell on 
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Satellite 



Figure 3.4-1. Rain cell division on a 
satellite communications 
link. 
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the electric field is determined by Eq. (3.3-4). The resulting fields 
after passing through the piecewise homogeneous rain then can be 
expressed by the multiplication of the rain cell matrices as follows: 





(3.4-1) 




Equation (3.4-1) is a very powerful tool in predicting the effects 
of an inhomogeneous rain on a millimeter wave communications link. 
However, data describing the true physical nature of a given rain 
extent are not readily available for most locations. To use this 
equation to its fullest extent. Sec. 3.6 presents a synthetic storm 
algorithm that can predict statistically rain effects for a variety of 
locations. Another use of this large scale discrete model is to allow 
only the first cell to contain ice particles. This is a physically 
meaningful situation because ice crystals usually exist only at the 
higher elevations of a storm. Physical situations of this nature are 
handled rather well by the above formulation and the generalized format 
of the scattering model has been achieved. 


3.5 A Discussion of Rain Reciprocity 


Now that a general model has been developed to predict the effects 
of an inhomogeneous rain on millimeter waves, the effects of nonreci- 
procity should be considered. Nonreciprocity as related to rain and 
millimeter wave propagation has received very little attention in the 
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literature*, however, under certain rain conditions it may have impor- 
tant consequences on millimeter wave communications links. To describe 
the effects of an inhomogeneous rain on millimeter waves, the medium 
is modeled by the product of M homogeneous rain cell matrices. This 
section will show that an inhomogeneous rain may be nonreciprocal 
and then the order of matrix multiplication in Eq. (3.4-1) is important. 

The reciprocity theorem [21] for electromagnetic fields sta^'es 
that the response of a medium to a source (transmitter) is unchanged 
when sdurce and measurer (receiver) are interchanged. Consider two 
arbitrary homogeneous rain cells that are placed adjacent to one another 
to form an inhomogeneous rain cell as shown in Fig. 3.5-1. The 
measured electric field after passing through the piecewise homogeneous 
rain medium can be expressed in matrix form as follows: 



The symbols A^, C^, and e. represent the entries for the rain 

cell matrix as defined by Eq, (3.3-1). For reciprocity as defined 
above to be satisfied, the matrices [S^] and [S^] must be cummutative. 
Using matrix algebra It can be shown that the two matrices are conmuta- 



1 

I 



Figure 3.5-1. An inhomogeneous rain, 
two cells of dissimilar 
canting angles. 
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tive and the medium reciprocal if and only if 


A . = 0 . 


(3.5-2) 


So that a more direct understanding of the nonreciprocal phenom- 
enon can be achieved, a simplifying assumption will be made. Assume 
that all the particles within the respective rain cells are equioriented 
oblate spheroidal raindrops with equal drop size (the distributions are 
represented by unit impulse functions). With this assumption, the 
matrix formulation of Eq. (3.3-1) is equivalent to the classical propa- 
gation constant model formulation [22]. The matrix entries and 

C- as defined by Uzunoglu, Evans, and Holt [5] are given below: 


A. = d|^ CDS^ 0^ + sin^ 0. 
'i ''i 

sin cos 8:. 


s- = 


h., 


sin^ 6- f d. 


cos-^ e. 


(3.5-3) 


where 




“jXA. N. fy 
e « « 1 


^h 


1 


"ju. N. fh- 
e 


(3.5-4) 


Using Eq. (3.5-3) and the definition of as defined by Eq. (3.5-1), 
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e. = C. - A. = (dy - ){cos2 0 ^. - sin^ e.} 

^ ^ (3.5-5) 

= Y^. cos^ Ze, . 

Now, Eq. (3.5-2) can be rewritten as 
y Y 

A = sin 2(02 - 9-]) (3.5-6) 

where 02 and 0 ^ represent the canting angles of the raindrops in the 
appropriate cells. The inhomogeneous rain medium of Fig. 3.5-1 is 
reciprocal only when A = 0 or when 

(02-9^)=^ , n = 0, 1,2,3... . (3.5-7) 

Note that the conditions under which reciprocity fails are independent 
of the length and the rain rate of the two cells, but only depend on 
the difference of the respective canting angles. 

By definition the classical propagation constant model requires 
that all drops in a giver thin, finite slab be equioriented. As a 
result, Eq. (3.5-7) indicates that the propagation constant model 
cannot deterministically model an arbitrary canting angle distribution 
without mathematically violating reciprocity. For this reason the 
3c2att&xn.ng model accounts for the canting angle distribution within 
the finite slab with the use of the generalized single particle 
scattering coefficients discussed in Sec. 3.1. 

Now the impact of nonreciprocity on communications link perfor- 
mance will be investigated. Although the assumptions of the classical 
propagation constant model lead to an exaggerated physical interpreta- 
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tion of the nonreciprocal phenomena, simple equations which explain the 
trends of nonreciprocity can be formulated. As a simple example consid- 
er the piecewise homogeneous medium of Fig. 3.5-1 with the incident 

* 4 

field being vertically polarized (E = 0). The exiting fields are found 

A 

• T 

from Eqs. (3.5-1) and (3.5-3) with Case 1 having the incident field E^ 
incident on Cell 1 and Case 2 having E' Incident on Cell 2. The fields 
for the two cases are given below: 


Case 1 

e(’* 2) = £A,B2 > H. e,B2)E; 

= (B^Bg + (A, + e^MAg + 52)}E^ 

Case 2 (3.5-8) 

i<2.1) . {A^Bj + B,A2 + ejBlJEy 

= {B^Bj + (A, + c,)(A 2 + e2>>Ey • 


Clearly 'if e-iBg equals e 2 B-| (a = 0), the fields of Case 1 equal the 
fields of Case 2 and the medium is reciprocal. 

Define aA to be the difference in attenuation of the incident 
field for the two cases. The general expression for aA with any inci- 
dent polarization is 


aA = 20 log 


10 


|0.2) 

?(2Trr 


dB 


(3.5-9) 


For a vertically polarized incident field, aA is 
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aA = 20 log-jQ 


i ( 1 . 2 ) 

fwj) 

y 


= 0 


(3.5-10) 


Independent of the rain conditions in either cell, attenuation is recip- 
rocal. One can show, using Eq. (3.5-9), that attenuation is recipro- 
cal for any arbitrary linear or circularly polarized incident field. 
Although it has not been proven for the general case of elliptical 
polarization, the same result is expected. Since attenuation should 
depend only on the volume of water present along the propagation path, 
the above conclusion seems reasonable. 

Now define AI to be the difference in channel isolation on a dual 
polarized communications systaa for the two cases of Fig. 3.5-1. The 
general expression for any arbitrary incident polarization is 


Al = 20 log.jQ 


£( 1 , 2 ) . £ 1 * 

|12.1) . !■* 


dB 


(3.5-11) 


-t- 

where E' represents the polarization state orthogonal to the incident 
field. For a vertically polarized incident field E’ = E x and Eq. 

A 

(3.5-11) reduces to 


Using Eq. 


Al = 20 log.|Q 
(3.5-8) in Eq. 


r(n2) 


(3.5-12) yields 


(3.5-12) 


Al = 20 log.^Q 


1 + 


A 

{B-jA2 + A.|B2 + ^ 


dB 


(3.5-13) 


Examination of Eq. (3.5-13) yields only two conditions under which 
reciprocity holds for isolation: 
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1) A « 0 

{B^Ag + A ^82 + 


The first case is when the entire medium is reciprocal* The second is a 
special case and will be Investigated shortly. 

In general isolation is nonreciprocal. This is to be expected be- 
cause it is the difference in the respective canting angles that makes 
the rain medium nonreciprocal. The magnitude of nonreciprocity depends 
on the difference of canting angles, input polarization, and other path 
rain conditions. The interdependence of the rain parameters {rain rate, 
rain extent, etc.} and their effect on the magnitude of Al will be 
presented through two special cases using the classical propagation 
constant model formulation. From these examples it will be evident that 
the magnitude of Al increases with frequency for a given input polariza- 
tion and path rain conditions. 

For the first special case allow the medium of Fig. 3.5-2 to have 
arbitrary rain rates in either cell. Allow the cells to have any 
length and JI 2 . The canting angle Sg is also arbitrary; however, let 
6^ equal zero. Under these conditions and using Eqs. (3.5-3) and 
(3.5-4) in Eq. (3.5-13) yield 


Al = 20 log-iQ 


N^) • ImCfy^ - 


dB 


(3.5-14) 


where 


0 -| - 0 
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, Cell 1 

Case 1 





Cell 2 


Figure 3.5-2a. Forward propagation. 



Cell 2 


Cell 1 

Case 2 
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Figure 3.5-2b, Reverse propagation. 


Figure 3,5-2, An example of non red procity. 
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02 t nu/2 , n = 0,1 ,2,3 .. . 

RR 2 0 
&2 0 • 

As long as 02 Is not zero or ninety degrees (or a multiple thereof), 
the magnitude of nonreciprocity is independent of the rain conditions 
in the second cell. At first this result may seem surprising; however, 
it is explained easily with a physical example. Figure 3.5-2a illus- 
trates Case 1. The incident field is not depolarized by Cell 1 and the 
drops of Cell 2 are oriented in such a manner that the field exiting 
Cell 2 is depolarized or rotated by <(i degrees. Case 2 is illustrated 
in Fig. 3.5-2b. In this case the incident field is depolarized or 
rotated by ^ degrees. The field incident upon Cell 1 is now oriented 
such that it will be depolarized. The depolarization by Cell 2 is the 
same in both cases since the input polarization is the same. However, 
the depolarizing properties of cell one are different for the two 
cases. Thus, the difference in isolation for the two cases depends 
only on the rain conditions of Cell 1. Equation (3.5-14) is plotted 
in Fig. 3.5-3 for three different frequencies as a function of rain 
rate. 

The propagation constant model formulation predicts serious non- 
reciprocal crosstalk effects on a given communications link using 
frequency reuse techniques. Although it provides simple equations to 
predict the effects of various rain parameters of rain reciprocity. 
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Figure 3.5-3, 

Theoretical predictions of the 
classical propagation constant 
model of the difference in 
isolation experienced on an 
uplink and a downlink of a 
millimeter wave satellite 
communications link. 
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the magnitude of the predicted results are exaggerated. The scattering 
model represents the rain in a more physically significant manner. With 
the inclusion of the spherical scatterer and the orientation and drop 
size distributions, the magnitude of aI is significantly decreased. 

Even for an operating frequency of 30 GHz, the scattering model indi- 
cates that AI is less than 0.5 dB. (The propagation constant model can 
also represent a drop size distribution.) 

The second special case will illustrate the effect that the 
difference in the respective canting angles (Sg - e-j) has on the 

magnitude of aI. For this example, the rain rates and path lengths of 

■>1 

both rain cells are equal, but the canting angles are allowed to be 
arbitrary. Under these assumptions equals and Eq. (3.5-13) 
reduces to 

AI = 20 log^Q 

where 

-j\N«,(f^ - f^) 

g - e 

for both cells. Note for the special case of 02 = ir - 9-| , Al equals 
zero. For any arbitrary linearly polarized incident field there is a 
certain special case of canting angle orientation that makes the medium 
reciprocal for isolation. This effect for vertical polarization is 
illustrated in Fig. 3.5-4. For a fixed value of 9-j , the isolation for 
each case and Al are plotted as a function of 02* whenever 

(Sg - 0^) " nir/2 where n = 0,1,2 Al is zero. Again, the magnitude 


1 + g tan a-j cot Ggl 
g + tan e-j cot e^ 


(3.5-15) 
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Figure 3,5-4. An example of the effect 

of canting angle on isolation 
nonreciprocity. (Classical 
propagation constant model , 
f * 3C GHz) 
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of AI is exaggerated. However, the peaks in Al for the special values 
of suggest that certain physical conditions may exist that could 
affect the operation of an adaptive isolation correction system on a 
given satellite communications link. 

The purpose of this section was to show that nonreciprocal condi- 
tions may exist on an uplink and a downlink of a given satellite commu- 
nications system. The interdependence of the parameters attributing to 
the nonreciprocal phenomena have been investigated with the use of 
simple examples. In general the scattering model does not seem to 
predict any significant nonreciprocal effects for typical rain situa- 
tions in the 10 to 30 GHz frequency band. Thus, the prediction of 
isolation nonreciprocity is exaggerated by models that do not include a 
canting angle or particle shape distribution. However, there still may 
be certain physical conditions that could have important consequences 
on millimeter wave communications links. Further investigation into 
the problem of rain reciprocity is necessary. 

3.6 Model Implementation 

Now that the scattering model has been presented formally and a 
word of caution has been noted about recip' ocity, the practical imple- 
mentation of the scattering model will be presented. The final goal of 
the scattering model is to predict signal attenuation, isolation, and 
phase shift for a variety of frequencies and site locations as a func- 
tion of measured rain rate data. Before this goal can be achieved, a 
detailed analysis of particle shape, particle size, particle density. 
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and particle orientation as related to ground rain rate must be 
formulated. These topics will be discussed in this section. Also, the 
relationship of the single particle scattering coefficients to eleva- 
tion angle S will be presented. The effect of non-ideal antennas on 
the predicted results will be investigated. Finally, this section 
presents a Rain Fropagation Prediation program (RFP) that includes a 
synthstio storm algorithm that can predict the effects of rain on milli- 
meter wave communications links for a range of frequencies, elevation 
angles, and locations. 

It has been shown that spherical raindrops, oblate spheroidal 
raindrops, and ice crystals may be used to model an arbitrary rain path. 
In tiie practical application of the tools developed in Sec. 3.4, the 
assumption will be made that a given cell consists entirely of ice 
crystals or entirely of raindrops. Tho first cell will usually be 
modeled by ice crystals and the remaining cells will be modeled using 
an effective percentage of oblate spheroidal and spherical raindrops. 

The distributions associated with the rain cell and ice cell will be 
presented separately, however, the assumption is made that both 
particle density distributions are a function of ground rain rate. It 
should be noted that the scattering model can include a mixture of rain 
and ice in a given rain cell, but presently data are not available to 
describe accurately this physical distribution of particles. The 
scattering model also will allow varying distributions along the propa- 
gation path. However, based on existing physica. data, each rain cell 
is assumed to have the same drop shape distribution, drop size distri- 
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bution, and drop orientation distribution. The drop density distribu 
tion is arbitrary for a given cell and is a function of rain rate, 

3.6.1 Raindrop Shape Distribution 

First consider the distribution of drop shapes. Jones [12] has 
found that the following distribution of shapes occur in an average 
rain for all radii greater than 0.9 mm: 


Spherical 

32% 

Oblate Spheroidal 

28% 

Prolate Spheroidal 

18.5% 

Irregular 

21,5% 


Since the scattering model is limited by the availability of single 
particle scattering coefficients, the above distribution will be 
modeled by an effective percentage of oblate and spherical raindrops. 
About one-half of the prolate drops would be expected to be aligned 
with their elliptical cross section facing the incident field and 
therefore produce depolarization. On a satellite link, the elevation 
angles are such that a fraction of the remaining prolate drops would 
also depolarize the incident field. The prolate drop class may be 
modeled by assuming the class is made up of 70% oblate drops and 30% 
spherical drops. A large fraction of the irregular drop class would 
also depolarize the incident field. If the irregular drop class 
is modeled by 90% oblate drops and 10% spherical drops, the total 
effective percentage of oblate drops is approximately 60%. Based on 
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the above arguments and good correlation with measured data (see 
Sec. 4.2) a 2 >cdn cell is assumed uo have 60<6 oblate spheroidal raindrops 
and 40X spherical raindrops. This distribution is assumed not to be a 
function of rain rate. 


3.6.2 Raindrop Size Distribution 


For a given rain rate there will be a distribution of drop sizes. 
There are several drop size distributions published in the literature 
[8], [23], [24], [25], A modification of the Laws and Parsons distri- 
bution will be used. The Laws and Parsons distribution as published is 
somewhat inconvenient in that it has no explicit functional form. 
However, real rains tend to deviate from any of the published distri- 
butions and an approximate functional representation of the Laws and 
Parsons distribution can be developed with little additional effect 
on the results. A simple approximation is that of a triangle. In 
Fig. 3.6-1, the Laws and Parsons distribution for a few rain rates is 
shown together with the triangular approximatirn where 


a 


a < 






0 el sewhere 

a fraction of drops in radius interval a* to a + da . 


V 


The symbol a^^ represents the modal or most frequently occurring drop 


Size distribution s(il 


- 47 - 



Figure 3.6-1. Triangular approximation to the 

taws and Parsons size distribution. 
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radios and a is the equi volumetric drop radius. Both radii are in 
inillimeters. The modal drop radius is a function of rain rate and is 
expressed by the relation [11] 

aj^ « 0.5 + 0.45 log.|Q RR (3.6-2) 

where RR is rain rate in mm/hr. The distribution of drop sizes can be 
represented by the relation 


n(a) 


Ny s(a) 


(3.6-3) 


number of drops per of space in the class a 
to a + d? . 


The symbol represents the particle density or the total number of 
drops per cubic meter of space. It can be shown (see Appendix 6.1.3) 
that Ny is a function of the modal drop size and thus is a function of 
rain rate, 

■ 5,833 RR a^^^ . (3.6-4) 

The calculation of the effective scattering coefficients (Eq. 
3.2-13) include summations over the single drop scattering coefficients 
weighted by the number of drops in the appropriate finite drop size 
interval. Single drop scattering coefficients are available typically 
at 0,25 mm drop radius increments up to 3.25 nin or 3.50 mm, depending 
on the frequency. By using curve fitting routines, the scattering 
coefficients can be determined at any point over the Interval 0.25 to 
3.50 nm. This permits the selection of any convenient drop size 


interval . 


In general define 
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» nwber of drop sise intervals over the range of drop 
sizes encountered 

a*^ ■ the equi volumetric drop radius value in m at the 

th 

midpoint of the 1 drop size interval 
Aa ■ drop radius interval in nw , 

Choose * 28 intervals. This provides a very fine division and the 
corresponding drop radius interval is aa ■ 0,125 m. The midpoints of 
these Intervals are 


a^ • 0.0625 + 0.125(1 - 1) 1 < i < 28 , (3.6-5} 

The number of drops per unit volume of space in the 1^^ drop size 
interval is 


a + ^ 

*1 *T 

NydJ « / . n(a) 

~ aa 

“ T 


(3.6-6) 


Using Eqs, (3.6-1) and (3.6-3) in Eq, (3.6-6) yield 

0.25 < a < djjj 


P" ®i 
m 


Ny(a^) • < 


0 


% i < 2«„ (3.6-?) 


«1 ‘ % 


A rough estimate of usable rain rates for this drop size distri- 
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bution can be formulated. For a less than 0.4 iwn the lack of small 

m 

drop size contributions would become significant. From Eq. (3.6-2) 
this modal drop radius corresponds to a rain rate of 0.60 mn/hr. For 
high rain rates the large drop contributions become more important. 

Drop radii out to 2a^^ are included in the triangular distribution. 

Since scattering coefficients are not complete beyond J = 3.25 mm, 

” 3.25/2 = 1.625 mm. From Eq. (3.6-2) this corresponds to a rain 
rate of 316 mm/hr. Thus, the above drop size distribution is reliable 
for rain rates between 1 mm/hr and 300 mm/hr. 

3.6.3 Raindrop Orientation Distribution 

The drop orientation distribution is very important to the predic- 
tion of cross-polarization. Although Brussaard [19] has shown that 
wind gradients affect the orientation of falling raindrops, little is 
known about the nature of the distribution itself. The work of 
Saunders [26] indicates that the distribution is Gaussian; however, 
researchers disagree on the range of the mean and the standard devia- 
tion. Most hav3 assumed an effective drop orientation and modeled the 
distribution with a unit impulse function. Although the unit impulse 
assumption leads to simple fcrmulations, it does not adequately 
describe the true physical situation. Uzunoglu, Evans, and Holt [5], 
using the distribution of Saunders, anploy a stochastic method of 
averaging the effects of a canting angle distribution; however, this 
is not internal to the model and does not consider the effects of arbi- 
trary orientation in the rain cell. Further experimental investiga- 
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tion into the problem of canting angle distributions is needed; 
however, until this is accomplished some sort of assumption must be 
formulated. Experimental measurements of cross-polarization at three 
frequencies over approximately the same propagation path (see Sec. 4.2) 
indicate that there is good correlation between theory and measurement 
with a Gaussian distribution whose mean ranges from -10“ to 10“ and 
with a standard deviation of 12“ . 

The drop orientation distribution within a finite slab can be 
discretized in the same manner as the drop size distribution. Define 

= number of drop orientation intervals over the canting 
angles encountered 

- the canting angle in radians at the midpoint of the 
orientation interval 

66* « canting angle interval in radians 

- mean canting angle in radians 

a = standard deviation in radians . 

Saunders indicates that canting angles ranging from -100“ to 100“ can 
be found in an average rain. Over this range choose Mj^ to be 20 
intervals. Thus, aJ is t/ 18 or 10 degrees. The midpoints of these 
intervals are given by the following equation: 

* -1.6581 + 0.1745 k radians (1 < k < 20) . (3.6-8) 


■f 
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th 

The number of drops per unit volume of space in the k orientation 
interval is 

N„ -{ — 

= — i- ; e " de (3.6-9) 

9.| 


where 


e, = le,^ - f} 
02 = te, + f } 


To evaluate the effective scattering coefficients of Eq. {3.2-13), 
the joint distribution of drop size and drop orientation must be known. 
Since the assumption has been made that the drop orientation is 
independent of drop size, the joint distribution is simply the product 
of the two distributions. 


• (3.6-10) 

3,6.4 Single Particle Scattering Coefficients and Elevation Angle 

It was indicated in Sec. 3,1 that the single particle scattering 
coefficients are a function of particle size, frequency and elevation 
angle. Coefficients are available for many frequencies and sizes; 
however, the selection of coefficients over the spectrum of e’’evation 
angles is limited. Uzunoglu, Evans, and Holt [5] have sho'-.r that a 
simple trigonometric relation can be used to evaluate the scattering coef- 
ficients for any arbitrary elevation angle. The relation is given below: 
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f (B) = fSPH g + f (0°) cos^ 6 

V V 

f^(6) = Sln2 S + f|j(0'’) cos^ e (3.6-11) 

where scattering coefficient evaluated for an elevation 

angle of zero degrees. This approximation is accurate to within 1% for 
elevation angles ranging from 0° to 90° and provides a simple formula- 
tion to evaluate scattering coefficients for a range of elevation 
angles. 

3.6.5 Ice Crystal Distribution 

Ice effects on satellite communications links have drawn consider- 
able attention since experimenters recently have observed cross- 
polarization without a significant fade in the co-polarized signal [27], 
[28], [29]. This so-called "anomalous" depolarization is attributed to 
a collection of ice crystals at high altitudes. Although the effect 
is not as severe, ice crystals also affect the isolation of a given 
communications system when followed by a series of rain cells. 

Haworth, McEwan, and Watson [30] have shown that the scattering 
coefficients for ice particles can be derived using Rayleigh scattering 
for frequencies up to 30 GHz, The derivation depends on the assumption 
that the particle eccentricity is very close to one. This physically 
models the ice particles as ice plates {oblate spheroidal) and ice 
needles (prolate spheroidal) as shown in Fig. 3,6-2 (see Appendix 



/ 

/ 


Ice Needle 


Figure 3.6-2. Ice particle shapes. 
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6.1.4). One disadvantage of this formulation is that it does not lend 
itself to a distribution of particle orientation. The assumption is 
made that all particles are equioriented at an effective canting angle. 
It has been observed by radar that lightning [31] and wind shear [32] 
align the particles in a given orientation. Since the ice cell is 
relatively thin, the assumption of an effective canting angle is not 
necessarily restrictive. The Haworth - McEwan - Watson formulation does 
include a particle size distribution. 

Using the Haworth - McEwan - Watson formulation in Eq. {3.2-13) 
along with the elevation angle formulation previously mentioned, the 
effective scattering coefficients of a thin finite slab of ice particles 
are given belov:.; 

^xx ° 'J' ^ + aJ®*-') + co$2 

S^y = -j ^ (cos2 9(A™ + A°®‘-') + sinA o(a[*'® + A^®*-')} (3.6-12) 

where 


aJ®*-' * sin2 s + cos2 3 A^®^ 
aJ®^' = sin2 3 A^*"” + cos2 s A^®*" 


(3.6-13) 


and A^^^, A^®^ and are complex constants independent 
of frequency and are given in Table 3.6-1 (see Appendix 6.1-4). Note 
that the scattering coefficients describing the prolate particle are 
independent of the elevation angle 6. The symbol V represents the 
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Table 3.6-1 » Ice scattering coefficients. 

j 0.026842} 
j 0.006179} 
j 0.026842} 
j 0.002674} 
j 0.009044} 


A™ = {6.812211 - 
= {3.268519 - 
= {6.812211 - 
aJ®*- = {2.150047 - 
A^*"^ = {3.95417 - 
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total volume of ice in a cubic meter of air. In Eq. (3.6-12) the 
assumption has been made that the ice slab has 50% ice needles and 50% 
ice plates. 

Haworth, McEwan, and Watson [30] use a value of V = 10“® in their 
predictions of anomalous depolarization. However, when a ice cell is 
followed by a series of rain cells, the volume of ice should have some 
dependence on ground rainfall. The hypothes-Cs is made that ike volume 
of iae per oubia meter of air in the freezing layer of a storm is equal 
to the volume of rain per auhia meter of air near the ground. The total 
volume of ice in a cubic meter of air is then a function of rain rate 
and is expressed by the relation 

V = (24,43 X 10-^) — . (3.6-14) 

m 

One should note that there still can be ice particles present without 
ground rainfall thus producing anomalous depolarization. 

3.6.6 Antenna Effects 

The development of the scattering model thus far has not included 
the interaction of the wave exiting the rain medium with the receive 
antenna. Even in clear weather conditions, a polarization mismatch is 
observed with non-ideal antennas. This polarization mismatch is 
modeled easily and for low rain rates has significartt effects in the 
prediction of attenuation, isolation, and phase shift. 

The co-polarized (main) and the cross-polarized (orthogonal) 
antenna polarization states can be described by a complex vector repre- 


sentation. The symbol e, represents the respective antenna state and 

a 

is defined below: 


X + y 


where for the co-polarized state 


(3.6-15) 


4 “ = costr^^} 


and for the cross-polarized state 


(3.6-16) 


= COS{Y } 

X cross'^ 


= sin{y }e 
y cross 


j 5 


cross 


The symbols y and 6 define the polarization ellipses of the respective 
antenna states as shown in Fig. 3.6-3. 

The wave-antenna interaction can be modeled by the use of a 
phasor voltage. Define a phasor voltage V(w, a) which describes the 
interaction of a wave polarization state w and an antenna polarization 
state as 


V(w, a) = ' V - (3.6-17) 

The wave state exiting the rain medium is found from Eq. (3.4-1) and 
the antenna states from Eq. (3,6-15). 

To evaluate attenuation, isolation, and phase referenced to the 
antenna ports, the clear weather wave state must be known, Tfie attenu- 



tan“^ { 


tan 

s1n 



{cos 2e cos 2t} 


cot"^ {axial ratio} 
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Direction of Propagation 


Figure 3.6-3. Polarization ellipse. 
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at ion of the co-polarized component is 


A = 20 log 


10 


V(w, 




= 20 Tog 


10 


'“w 

-> 

E ' 
w 


:co* 


CO’ 


dB 


(3.6-18) 


->• ■+• 

where is the clear weather wave and E^' is the wave exiting the rain 
medium. The isolation between the two receive channels is 


I = 20 1 


“SlO 


V(w'. a^^) 


V(w* , a T 
^ * cross ^ 


= 20 log 


10 


, ^ 

' 


:co’' 


• (3.6-19) 


Also of interest is the phase difference between the phasor voltages of 
the two receive polarization states co and cross. Define 


* = Cross Phase - Co Phase 
= Z{I„' -/_{%■ . a"*} radians . 


(3.6-20) 


The terms attenuation, isolation, and phase appear frequently in the 
remaining sections of this text and unless otherwise specified, the 
definitions of Eqs. (3.6-18), (3.6-19), and (3.6-20) will be assumed. 


3.6.7 A Synthetic Storm Algorithm 


Rain is not usually uniform over the extent of a storm and for 
increasing rain rates rain cells of higher rain rate tend to decrease 
in size. However, most researchers assume uniform rains and account 
for this phenomenon with an effective rain extent referred to in the 
literature as effective path length. As rain rate increases, effective 
path length decreases. The effective path length formulation is based 
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on experimental measurement's and theoretical predictions for uniform 
rain conditions. The intention is to formulate an effective path 
length based on statistical data that can be used for a variety of 
locations to predict rain effects on millimeter wave propagation. There 
is one major disadvantage to this approach: effective path lengths are 

frequency dependent. Thus, experimental data is required for all 
frequencies of interest and this is overly restrictive. The frequency 
dependence of effective path length is investigated in this section and 
a fpt'..p.nsnjy ind^psyi.i&Kt sunthstij alyovithtt that models nonuni- 

form rain rates is presented. 

The actual or measured attenuation on a given propagation path 
with physical rain extent L can be expressed as 
L 

= / A,{RR(0, f) di (3.6-21) 

m Q i 


where A^{RR(i), fl is attenuation per meter and is a function of rain 
rate and frequency. Note that rain rate is also a function of length 
Effective path length in km is defined as 


A 

^gf^^AVG* ~ A TRR fT 
e AVG ^AVG^'^^AVG* 


(3.6-22) 


where ^ theoretically predicted attenuation per 

kilometer assuming a uniform rain rate. The average rain rate 

is based on rain rate m’*asurements taken over a long period of time. 

The measured attenuation is also indirectly a function of average rain 
rate. Measured attenuation and measured rain rate data are compared on 
an equal probability of occurrence basis over a long time base (see 
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Sec. 4.1). This removes the Instantaneous time dependence of the 
measurements. Note that if rain rate is not a function of length £, 

\ ** ^AV6^*^*^AVG* *“* effective path length would equal the 

physical rain extent 1. 

To investigate the frequency dependence of effective path length, 
consider the ratio of two effective path lengths for twr frequencies f^ 


and fg: 


LfRR fl A fRR fl*^ A-i{RR(fi.), f-i}dA 
e^ AVG* r AVG^ AVG* 2' 0 


LgtflR^yG* ^2^ ^AVG^*^*^AVG’ ^1^ 


/ A,{RR(JI), fM 
0 ‘ ^ (3.6-22) 




Ratio is the ratio of measured attenuation for the two frequencies. 
Ratio Rp is a theoretical prediction of the ratio of attenuation 
assuming uniform ^ain conditions. For effective path length to be 
independent of frequency, Rp must equal inverse. The effective path 
lengths for two frequencies (19.04 GHz and 28.56 GHz) are plotted a;, a 
function of statistical (average) rain ra*te in Fig. 3.6-4. Effective 
path length data for these curves were calculated from three months 
(July, August and September, 1977) measured attenuation and rain rate 
data. The data were collected at the VPI&SU earth station using the 
COMSTAR D2 satellite. It is obvious from Fig. 3,6-4 that effective 
path length is frequency dependent. Other researchers have also 
observed this frequency dependence [33]. This indicates that for 
higher rain rates a theoretical model should not assume uniform rain 
conditions, but model the nonuniform rain with a synthetic storm 
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Figure 3.6-4. Effective path lengths for the 
VPI&SU COMSTAR 19 and 28 systems 
based on measured data collected 
during the months of July, August, 
and September, 1977. 
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algorithm. 

The scattering model, because of its ability to represent a piece- 
wise homogeneous rain, can accomodate easily a sy.thetic storm 
algorithm. The storm algorithm will model rain rate as a function of 
position along the path and storm extent as a function of elevation 
angle. Based on the data of Fig. 3.6-4, two physical assumptions will 
be made: 

1} the rain is uniform for lov^ rain rates 

2) as rain rate increases, the rain becomes nonuniform. 

The first step in the development of a synthetic storm algorithm 
is to obtain a functional relationship between effective storm extent 
and elevation angle. Define effective storm height as 

= L (10 m/hr) sin 3 (3.6-23) 

where L_ is the effective path length calculated for a rain rate of 10 
7m/hr , and 3 is the elevation angle. Since for low rain rates the rain 
is nearly uniform, effective path length at 10 mm/hr is frequency 
independent (R^ = and equal to the physical rain extent L. As a 
result, Eq. (3.6-23) can be rewritten as 

Hg = L sin 8 km (3.6-24) 

or the physical rain extent is 


km (8^0). 


(3.6-25) 
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Note that effective storm height approaches zero as elevation angle 
approaches zero. The effective storm height formulation provides a 
limit to the storm extant both in height and ground extent. Using the 
data of various experimenters [17], [34], [35] at various locations, 
elevation angles, and frequencies, the functional relationship between 
elevation angle and effective storm extent can be obtained. Valid 
statistical attenuation data are very limited over the spectrum of 
frequencies and elevation angles; however, enough data are available to 
provide a first order approximation for an effective storm model. The 
available data at a rain rate of 10 mm/hr are used to calculate effec- 
tive storm heights and these are tabulated in Table 3.6-2. Using the 
data of Table 3.6-2, an effective storm model is presented in Fig. 
3.6-5. From this figure given an elevation angle, the physical rain 
extent L can be calculated where L is not a function of rain rate. 

To model the nonuniform rain rate along the physical rain extent, 
assume that the maximum rain rate occurs in the rain cell closest to 
the receive antenna or cell M of Fig, 3.4-1. This is a valid assump- 
tion over long periocs of time since it has been shown that statistical 

rain rate is independent of the measurement location [35], The rain 

th 

rates of all other rain cells will be a function of the M cell. This 
functional relationship is the storm algorithm. 

To obtain the functional relationship, another comparison to 
measured attenuation data is necessary. However, for this comparison 
only data for a single frequency and a single location are needed. 
Figure 3.6-6 is a plot of measured attenuation data at 28.56 GHz as a 
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Table 3.6-2, Effective storm heights obtained from 
measured attenuation and rain rate 
statistics from various locations. 


Data Point 

3 

f 

He 

Location 

1) 

54° 

30 GHz 

5.4 km 

University of Texas 

2) 

44° 

28 GHz 

5.7 km 

VPI&SU 

3) 

44° 

19 lariz 

5.6 km 

VPI&SU 

4) 

23° 

28 GHz 

4.4 km 

Comsat Labs 

5) 

23° 

19 GHz 

4.4 km 

Comsat Labs 
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Figure 3.6-5, Effective storm model based 
on data from Table 3.6-2. 
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Figure 3.6-6. Measured 28 GHz attenuation 
ccxnpared to the theoretical 
prediction of the scattering 
model using the synthetic 
storm algorithm. (28.56 GHz) 
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function of statistical rain rate. The data were taken over a three 
month period (July, August and September, 1977) at the VPI&SU earth 
station. To remove the Instantaneous time dependence, the curve is 
plotted using equal probability techniques. For an elevation angle of 
44°, using Fig. 3.6-2 along with Eq. (3.6-25) the rain extent L is 
8.2 km. 

Assume that the rain extent can be modeled by ten rain cells as 
shown in Fig. 3.6-7. The rain rate of cells nine and ten equals the 
statistical rain rate of Fig. 3.6-6. Rain cells one through eight are 
equal and are related to cell ten by the following relation: 


RR^O ’‘i 

RR, = 


RR^q mm/hr . 


(3.6-26) 


This relation models the nonuniform rain rate over the rain extent. 

The value of x^. can be obatined by trial and error until the attenua- 
tion corresponding to a rain rate of 40 mm/hr is predicted by the 
scattering model. The value of Xg by definition is zero. The 
value of x^_^ is found to be 

x^^Q = -.66 . (3,6-27) 

Although the value of was obtained from one set of attenuation 
data, it models the nonuniformity of a rain medium as rain rate 
inc»'eases and is completely general for a variety of frequencies, 
elevation angles, and locations. 

The theoretical prediction for 28 GHz attenuation over the entire 
spectrum of rain rate is shown in Fig. 3.6-6. Other comparisons to 
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measured data using the above formulation are presented in Chapter IV. 
Given an elevation angle, the synthetic storm algorithm models non- 
uniform rain rates and effective stom extent and reliably predicts 
rain effects on communications links regardless of site location, 

3,6,8 Rain Propagation Prediction Program 

The Rain Propagation Prediction program (RPP) uses the scattering 
model and the synthetic storm algorithm of the previous section to 
predict attenuation, isolation, and phase as a function of rain rate. 
The input format allows changes in system parameters so that the 
effects of rain on millimeter wave propagation as a function of 
frequency, location and elevation angle can be predicted. The program 
is written in FORTRAN and its operational format is explained in 
Appendix 6,2. The program listing is found in Appendix 6.2.3. The 
RPP program is computationally efficient and it accurately predicts 
measured experimental data. 




CHAPTER IV 


EXPERIMENTAL VERIFICATION 

In this chapter the theoretical predictions of the scattering 
model are compared to measured data obtained from three different site 
locations: (1) VPI&SU, Blacksburg, Virginia; {2} Comsat Labs, 

Clarksburg, Mary' and; and (3) University of Texas, Austin, Texas, The 
data obtained from these sites cover the existing frequency range of 
interest for millimeter wave communications links (11 to 30 GHz) and 
provide a good comparison of theoretical predictions to measured data 
for a variety of system parameters. Before the data comparison is 
made a system description of the VPI&SU experiment is presented with an 
emphasis on general data reduction techniques. Then measured attenua- 
tion, isolation, and phase data for three frequencies (11.7 GHz, 19.04 
GHz, and 28.56 GHz) collected by the VPI&SU experiment are compared 
with theory. The scattering model predict'* ons are then compared to 
attenuation data measured by Comsat Labs and the University of Texas to 
demonstrate the flexibility of the synthetic storm algorithm. The 
chapter concludes with an investigation of frequency scaling techniques 
and the frequency dependence of attenuation and isolation. 

4.1 The VPI&SU Experiment 

4.1.1 General System Description 

The Virginia Polytechnic Institute and State University (VPI&SU) 
millimeter wave propagation experiment is sponsored by NASA with 
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important additive support from the Defense Communications Agency and 
the U.S. Army Research Office. The experiment monitors continuous 
transmissions from the Cofimuni cations Technology Satellite (CTS) at 
11.7 GHz and the COMSTAR D2 satellite at 19.04 GHz* and 28.56 GHz. The 
VPI&SU experiment is designed to measure attenuation, isolation, and 
phase shift and correlate these quantities with rain rate and other 
weather data. The experiment has been described in considerable detail 
in previous publications [36], [37], [17] and will be reviewed here 
only briefly with an emphasis on data reduction techniques. 

A block diagram of the VPI&SU system 1s presented in Fig. 4.1-1. 

The experiment is controlled by a PDP-11/10 mini computer which monitors 
the experiment through a digital controller and analog to digital 
converter interface. The three different frequency systems are moni- 
tored simultaneously aru‘ the PDP 11/10 stores receiver data along with 
weather and time reference data on a storage disk for later transfer to 
an IBM/370 for data reduction. The CTS (11.7 GHz) satellite transmits 
a right hand circularly polarized wave and the CTS ground antenna 
receives both the co-polarized and cross-polarized signals. The 
COMSTAR receive antennas are both dual -linearly polarized. The 
COMSTAR satellite transmits vertical polarization at 28.56 GHz; however, 
the COMSTAR 19.04 GHz beacon switches between vertical and horizontal 
polarization at a 1 KHz rate. This feature was incorporated into the 
propagation experiment to provide information on the differential 
attenuation and differential phase properties of the rain medium. As 
a result of this switching polarization, there are four 
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receive channels (Vertical Co, Vertical Cross, Horizontal Co, and 
Horizontal Cross) for the 19 GHz system. The 19 GHz IF signal 
processor has a deswitching control unit that provides a switching 
signal to monitor the various polarization states. The 28 GHz system 
is also a dual -polarized system and monitors the co-polarized and 
cross-polarized signals. From these eight channels, attenuation, iso- 
lation, and phase data are recorded for the different frequencies and 
polarizations. The dynamic nge of each system and other system 
parameters are tabulated in Tc^le 4.1-1. The polarization parameters 
describing the polarization ellipses of the incoming waves and antenna 
states are tabulated in Table 4.1-2. 

4,1.2 Preliminary Data Processing 

As indicated in the previous section, on-line data acquisition is 
controlled by the POP 11/10 it ni computer. Each data channel has an 
identifying number that lo ui throughout the data reduction process. 
The POP 11/10 also monitors vcrijus status indica ors *-hat indicate 
the validity of the data and whether it should be st' on the POP 
11/10 storage .disk. Data are stored only when the stci:us indicators 
indicate valid data and when the specific data input has changed by 
a predetermined amount. Data are stored with a time flag and other 
Identifiers for later off-line processing. 

During periods of low experimental activity, accumulated experi- 
mental data can be transferred to the main IBM 370 computer for sub- 
sequent data reduction. Before each data transfer all data are 



Table 4,1-1. VPI&Sl) system parameters. 


System Parameter 

(II-tIhz} 

COMSTAR 
{19.04 GHz) 

COMSTAR 
(28.56 GHz) 

Antenna 



■ 

Reflector size: 

12‘ 

4' 

4’ 

Gain; 

50.9 dBi 

44.7 dBi 

44.7 dBi 

Polarization: 

RHCP, LHCP 

LV, LH 

LV, LH 

Look angle 




Azimuth: 

240° 

213° 

213° 

Elevation: 

33° 

44° 

44° 

RF Front Ends 




Co {Gain/Noise Temp): 

41.5 dB/962° K 

34.8 dB/1540° K 

37.5 dB/1716° K 

Cross {Gain/Noise Temp): 

79 dB/184° K 

32.8 dB/1450° K 

34.5 dB/1716° K 

IF Signal Processors 




Dynamic Range {Co/Cross} 

39 dB/39 dB 

30 dB/30 dB 

39 dB/39 dB 

Clear Weather Input Power: 

-78 dBm/- 75 dBm 

-90 dBm/- 125 dBm 

-86 dBm/- 121 dBm 

{Co/Cross} 




Clear Weather Isolation: 

37 dB 

35 dB 

35 dB 
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Table 4.1-2. VPI&SU polarization parameters. 


Polarization Parameters 

CTS 

(11.7 3Hz) 

COMSTAR 
(19.04 GHz) 

COMSTAR 
(19.04 GHz) 

COMSTAR 
(28.56 GHz) 

Co-Polarized Wave 





e: 

-44.34® 

0.82® 

0.82® 

1.02® 

t: 

0® 

52.5® 

142.5® 

52,5® 

Co-Polarized Antenna 





e: 

-43.68 

0® 

0® 

0® 

t: 

0® 

52.5® 

142.5® 

52.5® 

Cross-Polarized Antenna 





e: 

43.68 

0® 

0® 

0® 

t: 

90® 

142.5® 

52.5® 

142.5® 


Note: See Fig, 3.6-3 for definition of e and t. 
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checked for time sequence and any required calibration is performed at 
this time. The data are transferred from the PDP-11 /TO storage disk to 
magnetic tape and then read by the IBM 370 onto a system data base 
tape. The data base tape is a aormlete time history record of the 
VPI5SU experiment. A data exception tape is stored also in the IBM 370. 
This tape records time markers of any system failure that cannot be 
determined by the status indicators such as a complete power failure 
and any other special cases in which erroneous data may have been 
stored. For a particular data reduction time interval , these two 
tapes are merged to form a process file. This process file can be 
accessed through an IBM time sharing system that allows efficient data 
reduction over the specified time interval. An inter-active graphics 
system is also available to aid the data reduction process. 

4.1,3 Data Reduction 

The data reduction software [38] is capable of processing data in 
a wide range of data reduction formats. A few of these formats will be 
illustrated in this section. Figures 4.1 -2a through 4.1 -2e illustrate 
rain rate and signal behavior during a storm observed on August 9, 

1977. These figures are time history plots and illustrate the instan- 
taneous time dependence of signal level and rain rate. Note that the 
rain began to fall after significant fad-'ng had occurr-'J at all three 
frequencies. The storm was severe enough 025 mm/hr) to cause the 
19 GHz and 28 GHz systens to lose phase lock. Note that the peak fades 
occurred approximately at the same time and the severity of the fades 
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RAIN RATE 
AUGUST S. 1977 


I 



TIME IN MINUTES 

Figure 4.1-2a. Ground rainfall rate on August 9* 1977. 

(The rain gauge is located beside the 
receiving antennas.) 


11.7 CO 
AUGUST 9. 19T7 



0 18 30 S4 72 90 

TIME IN MINUTES 

Figure 4.1-2b. CTS 11.7 GHz co-polarized signal 
behavior on August 9, 1977, 
(Measured fade 15 dB) 




19V CO 

AUGUST 9. 1977 


- ao - 


i 



Figure 4.1 -2c. 19 GHz vertical co-polarized signal 

behavior on August 9, 1977. 
(Measured fade 32 dB) 


28 CO 

AUGUST 9. 1977 


a 



Figure 4.1 -2d. COMSTAR 28 GHz co-polarized signal 


behavior on August 9» 1977. 
(Measured fade 32 dB) 




0 

Figure 4 


18 36 54 72 90 

TIME IN MINUTES 

l-2e, 11.7 GHz cross-polarized signal 

behavior on August 9, 1977. 


Figure 4.1-2. VPI&SU system time histcry plots. 
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increased with freqv 'y. 

Figures 4.1 -3a ana 4.1 -3b are scatter plots which provide an 
insight into the instantaneous interdependence of attenuation and 
isolation at different frequencies. Figure 4.1 -3a is a plot of isola- 
tion verses attenuation at 11.7 GHz for the storm of August 9, 1977. 
Figure 4,1 -3b is a comparison of attenuation of the 19 GHz vertical and 
the 28 GHz COMSTAR channels during the storm of August 9, 1977. The 
attenuation ratio for these two frequencies is approximately 2 to l. 
Attenuation ratios will be discussed in more detail in Sec. 4.4, 

The data reduction process used in Figures 4.1-2 and 4.1-3 did not 
include any statistical techniques. However, for the prediction of 
rain effects on millimeter wave signals the removal of the instantane- 
ous time dependence from measured data is necessary. This is accom- 
plished by using the VPI&SU reduction software coupled with a statis- 
tical analysis system called SAS [39]. System engineers designing 
conmuni cations links operating in the millimeter wave band are inter- 
ested in the percent of time that attenuation, isolation, and rain rate 
will equal or exceed a given value. Such "exceedence" plots are 
presented in Fig. 4.1-4 and 4.1-5. Figure 4.1-4 is a plot of the 
attenuation statistics for July, August, and September, 1977. Figure 
4.1-5 is a plot of the rain rate statistics over the same time base. 

The theoretical prediction of the Rice-Holmberg rain rate model [40] is 
also plotted in Fig. 4.1-5 and agrees rather well with the measured 
data. The Rice-Holmberg model will be discussed in more detail in 

i 

i 

Sec. 4.2.' 
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Figure 4.1 -3a. 11.7 GHz isolation versus attenuation 

for the storm of August 9, 1977. 


28 ATTENUATION, 08 



Figure 4.1-3b. A comparison of attenuation on the 

19 GHz vertical and the 28 GHz COMSTAR 
channels during the storm of August 9, 
1977. 


Figure 4.1-3. VPI&SU system scatter plots. 
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Figure 4.1-4. Measured attenuation statistics for 
July, August, and September, 1977, 
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Figu.-e 4.1-5. Rain rate statistics for July, 
August, and September, 1977. 
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In Fig. 3.6-6 a plot of attenuation versus rain rate was presented 
using equal probability techniques. This plot was obtained by compar- 
ing attenuation data and rain rate data of Figs, 4.1-4 end 4.1-5 for a 
given percent of time. This process removes the instantaneous time 
dependence as seen in Fig, 4,1-2 and represents the data on an equal 
probability of occurrence basis. Figure 3.6-6 is duplicated in Fig. 
4.1-6 for comparison to Figs. 4.1-4 and 4.1-5. 

4.2 Comparison of Theory with Measured Data from the VPI&SU Experiment 

This section demonstrates the ability of the scattering model 
to predict the effects of rain on a dual-polarized satellite communi- 
cations link. Theoretical predictions of the scattering model will be 
compared to measured attenuation, isolation, and phase data 
collected by the VPI&SU earth station during the months of July, 

August, and September, 1977. Experimental results are presented 
using graphs which also contain theoretical predictions of the rain 
propagation prediction program (RPP). These graphs make the compari- 
son between theory and experiment evident at a glance. The frequency 
range of the VPI&SU data and the polarization diversity of the VPI&SU 
experiment provide an excellent opportunity to test the scattering 
model with a variety of system parameters. 

Before theoretical predictions can be formulated for the three 
systems of the VPI&SU earth station, the input parameters for the RPP 
program must be determined. The only input parameters needed are the 
operating frequency, the elevation angle, the effective physical rail 



HTTENUHTION <DB) 



ts 

IS 

esr 

isr 

car 

IS 

G9 

tsi 

ta 

CS 

ta 


rn 

m 


RRIN 

RHTE 

( HM/HH ) 


Figure 4.1-6. 


Measured 28 GHz attenuation 
versus rain rate. {Plotted 
using equal probability of 
occurrence techniques.) 
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extent, and the system polarization parameters. The polarization para- 
meters and the elevation angles for the three systems are given in 
Tables 4.1-2 and 4.1-1, respectively. Knowing the system elevation 
angle, the ohysical rain extent L can be determined using Fig. 3.6-5 and 
Eq. (3.6-25). The effective rain extents for the CTS and COMSTAR 
systems are 10 Li and 8.2 km, respectively. The RPP program is 
currently limited to the discrete frequency values of 11.0, 14.0, 20.0, 
and 30.0 GHz. In the following figures the frequencies given are those 
of the measurement frequency and the particular frequency used in the 
associated theoretical prediction is the closest available value. 

4.2.1 Attenuation 

Engineers designing satellite communications links are faced with 
the problem of determining rain fade margins to provide a reliable 
communications link during periods of inclement weather. As the opera- 
ting frequency increases, rain fading becomes more severe. As a 
result, fade margins have to increase or ground station site diversity 
must be employed. Either of these solutions is expensive. Data are 
being collected to help determine the best solution to rain fade 
problems; however, for satellite links this collection of data is also 
expensive. For this reason, the development of a reliable theoretical 
model is necessary. The purpose of this section is to demonstrate the 
effectiveness of the scattering model in predicting the attenuation 
experienced on a satellite communications link for a given frequency 
and ground rain rate. 


- 89 - 


T^;e relationship between attenuation and ground rain rate as 
meas»i!red by one tipping bucket rain gauge near the VPXSSU antennas is 
il lustrated In Fig. Measured data are presented using equal 

probability data reduction techniques. The theoretical predictions are 
those of the RPP program where the nonuniform rain rates present on a 
satellite link have been modeled using the syni^ettc atom atgovi-bhn 
presented in Sec. 3.6. As can be seen, the agreement between theory 
and experiment is very good over the entire rain rate spectrum for all 
three frequencies. The severe fading at the higher frequencies (19.04 
and 28.56 GHz) indicate that site diversity may be the only solution 
for communications links operating at these frequencies and requiring 
a high degree of reliability. 

Simple equations relating attenuation to rain rate can be written 
in the form 

A « a RR^ dB (4.2-1) 

where A is attenuation in dB and RR is rain rate in mm/hr. Using power 
curve regression techniques on both the measured and theoretical data 
presented in Fig. 4.2-1, the constants a and b of Eq. (4.2-1) can be 
determined and are given in Table 4.2-1. These values include the 
effects of the synthetic storm algorithm and are valid for rain rates 
up to 60 mm/hr. 

Ground station site selection is aided by the knowledge of the 
percent of time that a given value of attenuation will be exceeded 
during a particular time interval. Thus, the percent of time that the 
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Figure 4.2-1. Attenuation versus rain rate for 
July, August, and September, 1977. 
(11.7, 19.04. 28.56 GHz) 
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Table 4.2-1. Power curve fits to 

attenuation versus rain rate 
data presented In Fig. (4.2-1). 
(RR < 60 nini/hr) 


Frequency 

HBBBi 

j Heasured(VPISSU) 

a 

b 

y »2 

a 

b 

^2 

n GHz 

0.9646 


i.qo 

1.9446 

0.5195 

0.92 

19 GHz 

2.587S 

0.5598 

m 

2.327 

0.5615 


28 GHz 

6. 4779 

0.4569 

m 

5.8033 

0.4828 

rm 
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link performance is degraded beyond acceptable system specifications 
can be determined. Using the scattering model coupled with the Rice- 
Holmberg rain rate model [40]* attenuation statistics can be predicted 
that agree rather well with the measured attenuation statistics of the 
VPI&SU earth station. The Rice- Hoi mberg rain rate model is an empiri- 
cal fit to a very large data base of measured rain rate data. The 
model predicts the percent of time that a given rain rate was exceeded 
based on the total rain accumulation over the time period in question. 
The Rice-Holmberg equation is given below: 

%T{RR) * Aj^^{0.038,,^ g-0.03RR ^ q 2(| - + 1.86 

(4.2-2) 


where is the total rain accmnuUfcion in mm times 100 divided by the 
number of hours in the accumulation period, is the ratio of accumu- 
lation of "thunder storm rain" to the total accumulation of rain, and 
RR is rain rate in mm/hr. The total accumulation of rain measured by 
the VPI&SU earth station was 153 mm for the months of July, August, and 
September, 1977 (2208 hours). As a result, 


Aj^^ = 6.93 mm/hr . (4.2-3) 

Since the time period In question is the thunder storm season, 

« 0.6 . (4.2-4) 


The Rice-Holmberg prediction is compared to measured rain rate statis- 
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tics In Figc The United States Weather Service observer 

(located within 8 km of the VPI&SU rain gauge) measured an accumulation 
of 418 m for the same time period. This Indicates that over short 
statistical time periods (less than a year), the local terrain may bias 
the observed statistics, particularly In mountainous regions. Kowever. 
over long statistical time periods the design engineer should be able to 
use rain accumulation data from th^’ USWS and predict rain rate and 
attenuation statistics for a given location E41]. 

Attenuation statistics for the months of duly, August, and 
September, 1977 can be predicted using the theoretical power curve 
regressions In Table 4.2-1 and the Rlce-Holmberg rain rate equation. 
Given a value of attenuation, the associated value of rain rate can be 
determined using the Inverse of Eq. (4.2-1) given below: 

RR • e“ * inii/hr . (4.S-5) 

Using Eq. (4,2-5) In Sq. (4.2-2), the percent of time that a given 
value of attenuation Is exceeded can be determined. Figure 4.2-3 
compares the theoretical attenuation exceedence plots to those m^suned 
at the VPI&SU earth station. As can be seen for rain rates above 10 
mm/hr (.U), experiment and theory agree rather well for 
frequencies. . 

4.2.2 Isolation 

As the channel capacity of sfngie-pbTarited^^^s^ 
tions systems become saturatsd, dual -polarized ctmimunlcatlons system 
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Figure 4«2-2» Measured rain rate statistics for 
July, August, and September, 1977 
compared to the Rice and Holmberg 
rain rate model . 
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Figure 4.2-3. Measured attenuation statistics for 
July, August, and September, 1977 
compared to the theoretical predic- 
tions using the Rice and Hti»mberg rain 
rate model and the scattering model . 
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have been suggested [1] to Increase channel capacity without signifi- 
cantly increasing the system cost. However, for operating frequencies 
above 10 GHz the isolation between the dual -polarized channels is 
affected significantly by rain. In this section theoretical predictions 
of isolation as related to attenuation are compared to measured data 
taken by the VPI&SU earth station during the month of August, 1977, 
Theoretical predictions of isolation statistics for the months of 
July, August, and September, 1977 are also presented along with the 
effaces of ice on the isolation of a dual -polarized communications link. 

In the previous section measured data were presented using equal 
probability data reduction techniques. However, multiple frequency 
isolation statistics are not presently available because of data 
reduction difficulties. Since hard rains often cause the COMSTAR 19 
and 28 systems to lose phase lock, uncertainties have developed in how 
best to handle these periods when computing isolation statistics. This 
problem was overcome in the evaluation of attenuation statistics by 
frequency scaling from lower frequencies where the receiver remained 
phase locked. However, as of this writing no standard format has been 
decided upon to scale isolation. Frequency scaling is discussed in 
more detail In Sec. 4.4. Measured isolation data presented in this 
section were reduced using instantaneous data reduction techniques. 

For a given value of attenuation, all of the corresponding values of 
Isolation are averaged and a standard deviation computed. For a given 
value of attenuation, the scatter in the associated isolation values 
is significant and can be explained by examining the physical scatter- 
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ing properties of the rain medium. 

The scatter in measured isolaticn data can be accounted for by 
three physical properties of the rain medium: 1) the oscillation of 
the mean canting angle with wind gusts [42], 2} the presence or the 
absence of a freezing layer, and 3) the movement of inhomogeneous rain 
cells along the propagation path. The assumption is made that wind 
gusts can cause thw mean canting angle to range between -10° and 10®. 
However, for a prevailing wind direction, the mean canting angle is 
predominately negative or positive [42]. Attenuation is relatively 
insensitive to changing canting angle; however, isolation can vary un 
to 3 or 4 dB depending on the canting angle and rain Intensity. The 
presence of an ice layer during a rain storm has no significant effect 
on signal attenuation; however, channel isolation degrades with an 
"increase in ice particle concentration. The effects of canting angle 
oscillation and the presence of an ice layer on the isolation of a 
dual -polarized communications link can be predicted by the scattering 
model. Figures 4.2-4, 4.2-5, and 4.2-6 compare the predictions of the 
RPP program to measured isolation and attenuation data during the month 
of August, 1977. In each of the fore-mentioned figures there are four 
theoretical curves a, b, c, and d with the following mean canting angle 
and ice content: 
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Figure 4,2-4. Isolation versus attenuation 
for the VPI^SSU CTS system for 
August, 1977, (11.7 GHz) 
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Figure 4.2-5. Isolation versus attenuation 

for the VPiaSU COMSTAR 19 system 
for August, 1977. (19.04 GHz) 
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Figure 4.2-6. Isolation versiis attenuation for 
the VPI&SU COMSTAR 28 system for 
August, 1977. (28.S6 GHz) 
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The ice layer Is assumed to be 1 km thick and since the ice layer does 
not affect signal attenuation significantly, the effective physical 
vcdn extent L ronains the same. Theory and experiment agree rather 
well for all three fY'equencies. However, close examination of the 
three figures indicates that the scatter in isolation increases with 
frequency. This is not completely accounted for by the oscillation of 
the canting angles or the presence of an ice layer. Because the data 
presented in the three figures are well within the system dynamic 
range, the scatter should not be a result of receiver error. One 
possible explanation is the movement of inhomogeneous rain cells along 
the propagation path. As indicated in Sec. 3.5, isolation may be non- 
reciprocal under inhomogeneous rain conditions and the magnitude of 
this nonreciprocity increases with frequency. So for a given value of 
attenuation (attenuation is reciprocal), there could be a range of 
varying isolation values depending on the inhomogeneous rai: conditions. 

Mean isolation can be calculated from attenuation through equa- 
tions of the form 

I = U + V In (A) dB (4.2-6) 

m 

where I is the mean isolation and A is attenuation in dB. If reliable 
m 

values of U and V can be determined for arbitrary earth station 
locations, then Eq. (4.2-6) could be used to calculate isolation 
statistics from existing attenuation statistics. This would enable 
communications engineers to predict the performance that would be 
available if an existing single-polarized link were replaced by a 
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dual "polaH zed link. The theoretical prediction assuming no Ice and a 
canting angle of -10** appears to best describe mean isolation as a 
function of attenuation for all three frequencies. With this assump- 
tion and using logarithmic curve regression techniques, the constants 
U and V of Eq. (4.2-6) can be determined for the mean Isolation data 
presented In Figs. 4.2-4, 4.2-5, and 4.2-6 and are given In Table 4 2-2. 
The logarithmic curve fits for this case for all three frequencies 
(11.0, 20.0, and 30.0 GHz) are presented In Fig. 4.2-7 to illustrate the 
frequency dependence of Isolation and attenuation. For a given value of 
attenuation, channel isolation Improves with frequency. It should be 
noted that the 11 GHz system Is circularly polarized and as a result 
1*^olat1on Is slightly worse than if the system were linearly polarized. 

Isolation statistics can also be predicted without a prior know- 
ledge of attenuation statistics. Using the scattering model coupled 
with the Rice-Holmberg equation. Isolation statistics can be predicted 
for arbitrary site locatiW based on total rain accumulation. Mean 
Isolation can be calculated from ground rain rate through equations 
of the form 

» T + W In (RR) dB (4.2-7) 

where Ij^ is the mean isolation In dB and RR is rain rate in mm/hr. 

The constants T and W again can be determined using logarithmic curve 
fitting techniques and are given In Table 4.2-3. For the theoretital 
assumptions of Fig. 4.2-7 (no Ice, "»'-10®), given a value of tsoTa** 
tion the associated value of rain rate can be determined using the 
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Table 4.2-2, Logarithmic curve fits to mean isolation 
versus attenuation for August, 1977. (For 
theory, A < 40 dB; for measured CTS and 
COMSTAR 19, A < 15 dB; and for measured 
COMSTAR 28, A < 30 dB) 



0y - 10° 

Theory RPP , 

No Ice 

Measured (VPI&SU) 

Frequency 

U 

V 

B 

U 

n 

r2 

11 GHz 

42.2507 

-9.4005 

1.00 

44.74 

-9.8 

0.76 

19 GHz 

40.1418 

-7.6020 

1.00 

46.56 

-10.63 

0.7087 

28 GHz 

42.6301 

-7.0897 

0,99 

39.38 

-6.67 

0.2469* 


* r^ is low because of the scatte*' in the data. 
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Table 4.2-3. Logarithmic curve fits to theoretical 
isolation versus rain rate data, (Rain 
rate < 60 nm/hr, e - -10®, no ice) 


Frequency 

|[HH 

W 

|h2 

11 GHz 

42.7801 

-5.8524 

1.00 

20 GHz 

32.7645 

-^.,2008 

0.99 

30 GHz 

29.0694 

-3.1824 

0.98 
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inverse of Eq. (4.2-7) given below; 

RR = eJ(p{(I|n - T)/W} m/hr . (4.2-8) 

Using Eq. (4.2-8) in Eq. (4.2-2), the percent of time that isolation is 
less than a given value can be determined. Theoretical isolation 
statistics are presented in Fig. 4,2-8 without validation by experiment 
for all three frequencies. Isolation statistics are affected by para- 
meters other than rain such as tracking error and antenna polarization 
sensitivity. Figure 4.2-8 only accounts for rainfall degradation. 

A decrease in channel isolation also can be observed without an 
associated rain event. This phenomenon is due to ice crystals which 
can depolarize incident fields without significantly affecting the co- 
polarized signal. Assuming that a 1 km freezing layer can exist 
without ground rainfall. Fig. 4.2-9 illustrates the effect this 
freezing layer would have on the VPI&SU earth station for various 
particle canting angles. Although these effects have not been observed 
by the VPI&SU facility, ice depolarization has been observed in other 
climatic zones [29]. The theoretical data in Fig. 4.2-9 is presented 
to illustrate the ice depolarization phenomenon and to demonstrate 
the flexibility of the RPP program. Note that the isolation of the 
circularly polarized 11 GHz system is relatively constant over the ice 
crystal canting angle range. 

4.2.3 Phase 


To Improve the channel isolation of a dual -polarized communica 


* TIME I1BSCI59I LESS THHM OR EQIBLLEI> 
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Figure 4.2-8. Theoretical prediction of the 
Isolation statistics for the 
VPISSU system. (July, August, 
and September, 1977) 
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Figure 4,2-9. Theoretical predictions of the 

scattering model of ice depolari- 
zation in the absence of rain for 
different particle orientations. 
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tions link during precipitation events, static or adaptive cancellation 
systems are being considered. These systems adaptively change the 
polarization state of the receive antenna to Improve system isolation. 
To design adaptive cancellation systems, a prior knowledge of the rela- 
tive phase of the dual -polarized channels is often helpful. For 
convenience the term ’’phase”, represented by the symbol 4, is defined 
as the cross-polarized channel phase minus the co-polarized channel 
phase. By controlling the relative phase $ between the two channels 
and the differential channel attenuation, system isolation can be 
improved. This section illustrates the sensitivity of phase changes to 
changes in antenna polarization parameters and the ability of the 
scattering model to predict these changes as systan isolation degrades 
with an increase in ground rain rate. 

Theoretical phase data versus isolation data are ccmipared to 
measured data in Figs. 4.2-10 and 4.2-12. Thp theoretical predictions 
assume that an ice layer Is not present along the propagation path and 
that the mean canting angle ( 0 *^) ranges between -10* and 10* as a 
result of wind gusts. Since adaptive cancellation syst^s correct 
for relatively fast changes in phase, measured data are presented on 
an instantcmeous time basis for a particular rain storm. This elimi- 
nates slow phase variations over longer periods of time and more 
accurately describes the requirements of an isolation enhancement 
system. 

Figure 4.2-10 represents a typical response of the CTS {11.7 GHz) 
system during the storm of August 9, 1977. Note that the. scatter in 
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Figure 4.2-10. CIS isolation versus 
phase for the storm 
of August 9, 1977. 

(CIS is circularly polarized) 
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isolation and phase supports the assumption of a 20** range in canting 
angle. The boxed-in area of Fig. 4.2-10 represents the system behavior 
before the actual rain event. This enhancement of isolation before the 
onset of a storm is seen often and is attributed to the special 
characteristics of the system polarization parameters. The circular 
polarization of the CTS system also accounts for the relatively large 
phase scatter for a given value of isolation. For circular polariza- 
tion, Overstreet [43] has shown that phase scatter is usually twice 
that of the mean canting angle range and this is verified in Fig, 
4.2-10. 

Phase changes during precipitation events for a dual circularly 
polarized communications syst«n are very sensitive to changes in 
antenna polarization parameters. If the satellite and ground station 
antenna polarization parameters are changed slightly from those found 
in Table 4.1-2 (the new epsilon's and tau's are given in Table 4.2-4), 
the phase retardation of Fig. 4.2-10 can change direction and advance 
as seen in Fig. 4.2-11. Note that for high values of isolation or for 
low rain rates, the phase has changed by a significant amount (61®) 
relative to Fig. 4.2-10; however, for higher rain rates the change is 
relatively minor (7®). For high rain rates, changes in the polariza- 
tion parameters of the spacecraft antenna are masked because of the 
large cross-polarized component generated by the rain. Changes in the 
ground station antenna are then minor for small polarization parameter 
deviations. However, for low rain rates the Interaction of the space- 
craft and ground station antennas becomes significant. The co-polar- 
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Table 4.2-4. Epsilon's and Tau's used in Fig. 4.2-11, 
(Clear weather isolation is 38 dB) 


Polarization Parameters 

CTS 


-43.35 

■'w= 

0° 

^co* 

-44.01 

^co* 

0° 

^cross' 

-44.11 

^cross' 

go'* 
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Figure 4.2- 


n. Theoretical prediction 
of isolation versus phase 
for the polarization 
parameters of Table 4.2-4. 
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ized phase is relatively insensitive to changes in the antenna polariza 
tion parameters; however, for low rain rates the cross-polarized phase 
is very sensitive to polarization parameter deviations. During low 
rain rate conditions, the cross-polarized component generated by the 
satellite is approximately equivalent to the cross-polarized component 
generated by the rain. Thus, the phase of the resultant wave is very 
sensitive to the cross-polarized phase of the wave transmitted by the 
satellite. Because phase {$) is insensitive to antenna polarization 
changes during high rain rate conditions, phase for a circularly 
polarized system can advance or retard depending on the antenna polari- 
zation parameters and the associated effect on the cross-polarized 
phase for low rain rates. 

Although the same arguments apply for a dual linearly polarized 
comnuni cations system as for the circularly polarized antenna systan, 
the uniqueness of the polarization parameters for a linear system 
reduce the dependence of phase on polarization parameter deviations. 
Also, the associated phase scatter of a linearly polarized system is 
significantly reduced as seen in Fig. 4.2-12. This figure represents 
a typical response of the COMSTAR 28 (28.56 GHz) systan during a storm 
on August 30, 1977. From the data presented In this section, a 
linearly polarized systan is recommended over a circularly polarized 
systan when an adaptive polarization systan is being considered [43]. 

4.3 Comparison of Theory to Measured Data from Other Experiments 


In the previous section the theoretical predictions of the 
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Figure 4.2-12. COMSTAR 28 Isolation versus 
phase for the storm of 
August 30, 1977. (COMSTAR 
28 is linearly polarized) 
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scattering model were compared to measured data from the VPI&SLt earth 
station for a variety of system parameters. Excellent correlation 
between measurement and theory was seen. However, the true test of a 
theoretical propagation model Is the ability to predict th& effects of 
rain on dual»po1ar1zed comnunl cations links for different frequencies, 
elevation angles, and Zoaaticna* Using the first order approximation 
of the synthetic stonn model (Fig. 3.6-5) and the respective system 
parameters of the Comsat Labs COMSTAR experiment and the University of 
Texas ATS-6 experiment, theoretical attenuation predictions of the 
scattering model are presented 1n this section that agree rather well 
with published measured data from those experiments. 

Since attenuation is relatively insensitive to changes In polari- 
zation parameters and the exact polarization parameters for the 
experiments of Comsat Labs and the University of Texas are net known 
by the authors, the polarization parameters of the VPI&SU COMSTAR 
system will be assumed In the theoretical predictions of this section. 
Comsat Labs monitors continuous transmissions from the COMSTAR D1 
■satellite at an elevation angle of 23®. As a result, using Fig. 3.6-5 
and Eq. (3.6-25) the effective rain extent Is 11.2 km for the 19.04 GHz 
and 28.56 GHz COMSTAR D1 beacons. Putting this data into the RPP 
program yields the theoretical attenuation predictions in Fig. 4.3-1 
for the 19 GHz and 28 GHz systems of Comsat Labs. The associated 
measured data were taken from published rain rate and attenuation 
statistics for July, 1976 through January, 1977 [34]. As can be seen, 
the agreenent between theory and measurement Is very good. Since no 
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Figure 4.3-1. Attenuation versus rain rate, theoretical 
predictions of the scattering model 
compared to measured data of Comsat Labs. 
(19.04, 28.56 GHz) 
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Isolation data are available* no theoretical predictions are presented. 
The University of Texas monitored transmissions from the ATS-6 
satellite on a time-available basis from July, 1974, to Hay, 1975. 
Measured attenuation data from the NASA Technical Note [35] along with 
the theoretical predictions of the scattering model are presented In 
Fig. 4.3-2. The ATS-5 propagation frequency was 30 GHz and the 
elevation angle from the University of Texas experiment site was 54®. 

As a result, the effective rain extent was 6.7 km. Again there is good 
agreement between theory and measurement. 

Figures 4.3-1 and 4,3-2 demonstrate the flexibility of the synthe- 
tic storm algorithm of Sec. 3.6 and the scattering model In predicting 
rain effects on millimeter wave conmuni cations links. However, before 
the scattering model can be used to predict the effect of rain on my 
mbitraxy earth station, more data must be collected to improve the 
synthetic storm model of Fig, 3.6-5. An effective storm height of 
6 km and a ground rain extent of 10 ton seems reasonable; however, 
these values may change for different climatic zones. Before the full 
potential of the synthetic storm algorithm can be achieved, more data 
are needed, 

4.4 Frequency Scaling 
4.4.1 Attenuation 

It is often very useful to scale measured attenuation data at 
one frequency in order to estimate system performance at another 
frequency. The frequency dependence of attenuation and the relation- 
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Figure 4.3-2. Attenuation versus rain rate. 

theoretical predictions of the 
scattering model compared to 
measured data of the University 
of Texas. (30 GHz) 
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ship to ground rain rate 1s predicted easily by the scattering model . 

In this section theoretical predictions of the scattering model are 
compared to measured data and to the recent frequency scaling algorithm 
of Hodge [44]. 

Figures 4.4-1, 4.4-2, and 4.4-3 illustrate the frequency depen- 
dence of attenuation measured at the VPI&SU earth station and the 
corresponding predictions of the scattering model. Measured and 
theoretical data are plotted using equal probability of occurrence 
techniques. Figure 4.4-4 presents the measured attenuation data of 
Comsat Labs along with the associated prediction of the scattering 
model. Although the rain extent and the elevation angles are different 
for the two sites, the ratio of the 28 8Hz attenuation to the 19 GHz 
attenuation Is the same. 

For a given value of attenuation at a particular frequency, the 
corresponding scaled value of attenuation at a second frequency can be 
determined using equations of the form 

Ag-uA^i dB (4.4-1) 

where A^ represents the attenuation of the 1 frequency in d8. Using 
power regression techniques the constants u and v In Eq. (4.4-1) can 
be determined and are given In Table 4.4-1 for both the theoretical 
and measured data presented in Figs. 4.4-1, 4.4-2, 4.4-3, and 4.4-4, 

Another way of describing the frequency dependence of attenuation 
Is through equations of the form 
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Figure 4.4-1, 28 GH 2 attenuation versus 

19 6Ht attenuation (VPI&SU). 
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Figure 4.4-2. 19 i^z attenuation versus 

n GHz attenuation (VPUSU). 
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Figure 4.4-3. 28 GHz attenuation versus 

n GHz attenuation (VPI&SU). 
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Figure 4.4-4. 28 SHz attenuation versus 

19 GHz attenuation (Comsat Labs). 


Table 4.4-1. Power curve fits to attenuation versus 
attenuation data presented in Figs. 

4.4-1, 4.4-2, 4.4-3, and 4.4-4. (A. < 40 dB) 

See Eq. (4.4-1). ’ 
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^2 _ ^2 
A ^ KK 


(4.4-2) 


where RR is rain rate in mra/lir and and b. are found in Table 4.2-1 
for the VPI&SU experimental systen. It is obvious from Eq. (4.4-2) 
that the frequency dependence of attenuation is also dependent on 
ground rainfall. 

Hodge [44] also has considered the dependence of attenuation 
scaling on rain rate. Assuming that rainfall along the propagation 
path is a Gaussian function of position on the path» the ratio of 
attenuations at two different frequencies is expressed by the relation 
A, 


“2 _ /> .. r.rSh ■ ^ 1 ^ 


no Un f ' ■- V 


(4.M) 


where and are cor Vv ■ p^sociated with power curve fits to 
attenuation per kilomett ' . ific attenuation)’ as a function of a 

uniform rain for the i^" frequency. These constants are given in 
Table 4.4-2. The scattering model prediction, the theory of Hodge, 
and a curve fit to measured data for the ratio of 28 GHz attenuation 
to 19 GHz attenuation are illustrated in Fig. 4.4-5. The theory of 
the scattering model and the theory of Hodge both employ a nonuniform 
distribution of rain rate along the propagation path. As a result 
both models agree rather well with measured data. However, unlike the 
Hodge formulation, the synthetic storm algorithm is internal to the 
scattering modal and can be used also to scale isolation (the Hodge 
formulation can predict only the ratio of two attenuations). 



Table 4.4-2. Power curve fits to theoretical 
attenuation per kilometer 
versus rain rate assuming uniform 
rain conditions. 


Frequency 

tt 

C 

r^ 

n 

0.0217 

1.2001 

0.98 

19 

0.089 

1.094 

0.99 

28 

0.276 

0.903 

0.99 
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Figure 4.4-5. Theoretical attenuation scaling, 
comparison of the Hodge and 
the scattering model formulations 
to measured data. 



4.4.2 Isolation 


Frequenev scaling of isolation is more complicated than scaling 
attenuat'*- I. Isolation is very sensitive to changes in rain conditions 
and as a result significant scatter in isolation occurs. Isolation 
measurements are also very sensitive to antenna polarization parameters 
and tracking errors. As a consequence of the large variation in 
isolation for a single frequency* the scaling of isolation presents a 
difficult problem to communications engineers* researchers and experi- 
menters. 

For a particular set of system parameters, the scattering model 
can predict the frequency dependence of mean isolation. Figure 4.4-6 
illustrates the frequency dependence of mean isolation for the system 
parameters of the VPI&SU earth station. The data presented in Fig. 
4.4-6 can be represented with equations of the form 

L t I dB (4.4-4) 

m2 m-^ 

where is the mean isolation in dB for the i frequency. The 
symbols t and w are constants containing the effects of the particular 
system parameters (canting angle, elevation angle, rain extent, etc.) 
on mean isolation. These values are tabulated in Table 4.4-3 for the 
VPI&SU experimental system. Although the data In Fig. 4.4-6 are 
presented without experimental verification, an insight into Isolation 
scaling has been obtained. 
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Figure 4.4-6* Theoretical mean isolation 
versus mean isolation for 
11. 20. and 30 GHz. 



Table 4.4-3. Power curve fits to 

isolation versus isolation 
data presented in Fig. 4.4-6 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

The effects of rain on milUmeter wave propagation will influence 
the future design of satellite communications systems. Therefore* it 
is necessary to predict the impact of precipitation on satellite 
communications system performance prior to the specification of system 
parameters. To this end a theoretical model has been developed to 
predict the degradation in system performance due to precipitation. A 
new deterministic model has been presented that models an inhomogeneous 
rain by discretizing the rain into several piecewise homogeneous rain 
cells. The scattering properties of each homogeneous rain cell are 
described by effective scattering coefficients. These coefficients 
describe the changes in polarization, attenuation and phase shift that 
an incident field experiences as it propagates through an ensemble of 
particles with a distribution of particle type {rain or ice), particle 
shape, particle size, particle density, and particle orientation. 

A frequency independent synthetic storm algorithm also was developed 
to model the effects of nonuniform rain rates present on satellite 
communications links. As a result, the acattsTing modal can predict 
accurately the effects of rain on a millimeter wave communications 
link for a variety of freruencies, elevation angles and locations. 

After a detailed derivation of the scattering model, the predic- 
tions of the associated Bain Propagation Prediction program were 
compared to available measured data. The comparison of theory to 
experimental data covered a wide range of systan parameters and 
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various site locations. Excellent agreement was obtained. 

Although the scattering model has predicted accurately the 
observed rain effects on satellite link performance at various ground 
terminals 3 more experimental data are needed to further refine the 
synthetic storm algorithm. Data describing the distribution of canting 
angles and particle shape along the entire rain extent also would lead 
to further refinements of the physical modeling of the rain conditions 
present on a satellite link. After these refinements and further 
verification of the model with a more extensive experimental data 
base, the scattering model vrill provide the system design engineer 
with a low cost, reliable model to predict the effects of weather 
on arbitrary satellite communications systans. 


CHAPTER VI 


APPENDIX 


6.1 Derivations 

6.1.1 Derivation of the Generalized Single Particle Scattering 
Coefficients 


In the development of the scattering model, a generalized single 
particle scattering coefficient was used to allow particles within a 
thin rain slab to have any arbitrary orientation. Thus, a canting 
angle distribution within the rain slab is possible. The derivation 
of the generalized single particle scattering coefficients is 

straight-forward and based on simple geometric principles [11]. 

*i 

Consider Fig. 6.1-1. An incident field E aligned with the 

arbitrary x-axis is incident on an arbitrary oblate particle canted 

*i 

at an angle 9 with respect to the x-axis. E' is decomposed into 

A 

*i 'i 

components E^ and Ej^ along the principal axes (v, h) of the oblate 
particle; respectively, 


• i ‘“i 
E^ = e’ sin e 

eJ = eJ cos e . 


( 6 . 1 - 1 ) 


*i *i *s *s 

Ey and are scattered by the oblate particle producing EJ and E^, 


K ^ sin s 


Eh = eJ B 


( 6 . 1 - 2 ) 
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where and are the complex scattering coefficients of the particle 
minor and major axis, respectively (f^ and are defined in Sec* 3.1), 
eJ and have components and along the x-axis and components 
EJy and along the y-axis as a result of the rotation of the inci- 
dent field by the oblate particle and are given by 

^ « eJ sin 6 = sinS 0 

s " *’h » 

(6.1-3) 

COS e = f ^ eJ sin 6 cos a 

^hy * ^h ® * ^h 9 cos e 
The total scattered field component along the x-axis is 


E® - P + E? 
S ^vx *^hx 


« sin2 e + f cos^ 6} 


(6.1-4) 


The total scattered field component along the y-axis is 


E® « E® - E? 
y vy hy 

* E^(fy - fjj) sin e cos e . 

The scattering coefficient is defined as the ratio 
E^ 

f^j^ ® * fy sin2 0 + f cos2 e 

^x 

The scattering coefficient f is defined as the ratio 

xy 


(6.1-S) 
of kl to il 


of E® to eJ. 
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^y/*i "" e cos e . (6.1-6) 

The coefficients f and f are obtained by replacing e in Eqs. 
yx yy 

(6.1-5) and (6.1-6) with 0 , where 
6 ^ = 9 - 90° 

and then 

V ^ ^xy ^ ® ® (6.1-7) 

and 

fyy = fy cos2 0 + f sin2 0 . (6.1-8) 

Equations (6.1-5), (6.1-6), (6.1-7) and (6.1-8) define the generalized 
single particle scattering coefficients as given in Eq. (3,1-1). 

k 

6.1.2 Evaluation of 2 g 
^ 

Before the effective scattering coefficients of Eq. (3.2-13) can 

k 

be determined, the evaluation of the summation s *9^ is necessary. 
The summation' extends over all the particles within the finite rain 
slab of Fig. 3.2-1 in the i^*^ size interval and the orientation 
interval. If the particles in the i and k intervals are numerous, 

the summation can be expressed as an integral, 

N. AA 277 « 

' ’ vii ° i i i ^ 

where n^^j^ is the number of particles per unit volume in the size 
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interval and the orientation interval. Performing the integral 
yields 


N 


i.Jc 


A& 2ir « -1 -j k„ pV2z 


9|li 0 0 0 ^ 


pdpd(^dz 

AA “* . -j k^ pS/22 

« 2tt n. ^ / 1 {p e ® } 


0 0^ 


dpdz 


A& ^ 


= _,• 2tt A& 
'o 


° ’J nrr^.k 


= -j m n. ' (6.1-10) 

where the following indefinite integral [45] was used: 

/ X e ^ ^ d 2 = ^ . 

Equation (6.1-10) is the result given in Eq. (3.2-7). 

6.1.3 Derivation of Ny 


Ny is the total number of raindrops per unit volume within a thin 
rain slab and is a function of rain rate and the distribution of 
drop size. Since drops of different size fall with different terminal 
velocities, they will contribute differently to the total measured 
ground rain rate. To account for this difference consider a time 
interval At over which all drops striking the ground will be summed to 
obtain the total rain accumulation. The drops in the drop size inter- 
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val a" to a + da have a terminal velocity of v^Ca) given below: 

v^(?) = 4.6 2a m/sec (6,1-11) 

where aT is in mm [46] . The height above the ground corresponding to 
the height of the last drops in the time interval At to strike the 
ground is 


hCa) = Vj(^At m 

= 6.5054VT At m (6.1-12) 

where H is in mm, Vj is in m/sec and h is in m. The volume of water 
accum' 'ated on top of the ground per unit area of ground surface in the 
time interval At due to drops in the size interval a* to if + da is then 

Vol, (a) = h(a) ^ (a)3 n(a)da (6.1-13) 

where n(a) is the number of drops per cubic meter of space in the size 
class ¥ to ? + d¥. This volume of water equals the rain rate for the 
drops with drop size ¥ to ¥ ■(■ d¥, RR(¥)d¥, times the time interval. 
Using Eq. (6.1-12) in Eq. (6.1-13) and eliminating At yields 

RR(¥)d¥ = v^(¥) ^ (¥)3 n(¥) d¥ , (6,1-14) 

The total rain rate is found by summing over all drop size contribu- 
tions, 

RR = ; RR(¥) d¥ 

= / v.p(¥) ^ (¥)3 n(¥> d¥ , 


(6.1-15) 
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Using the modified Laws and Parsons drop size distribution discussed 
in Sec. 3,6 gives 

RR « ^ (6.5054) Nw / (a)^^^ s(a) da 
^ ^ 0 

=* ^ (6.5054) Ny{4(2®/2 - 1)0/9 - l/lDa^'^^} 

* 47.624 Ny m/s (6.1-16) 

where is the modal drop radius in mm. Rain rate should be in units 
of mm/hr, but the above expression arose from terminal velocity in m/s, 
drop volume in mm^ and number density in per of space. Thus, to 
convert Eq. (6.1-16) into units of mm/hr, 

nr. _ /1.T » .7/2 0 0^ mm. r 3600 s,,10"3 m,, 

RR = 47.624 Nw a ' { — = — }{ — — }{ — 

V m m hr mm 

= 0.17144 Ny mm/hr . (6.i-17) 

To find the total number of raindrops per unit volume, the inverse 
of Eq. (3.6-17) yields 

Ny = 5.8327 (6.1-18) 

which is the result given in Eq. (3.6-4). 

6.1.4 Derivation of Single Ice Particle Scattering Coefficients 


In the derivation of the single ice particle scattering coeffi- 
cients, the theory of Rayleigh scattering can be used for frequencies 
up to 30 GHz [30]. Assume that an incident electric field is propaga- 
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ting in the 2 direction and is incident on an ice particle in the hv 
plane as shown in Fig, 6.1-2. The ice parti jle can be spherical, 
oblate spheroidal, or prolate spheroidal in shape. A dipole moment is 
induced in each of the principal planes of the particle. The assump- 
tion is made that the size of the particle is sufficiently small to 
ignore higher multioole moments. Each of the components of the inci- 
dent field along the principal planes {h, v) are considered to act 
separately in the excitation of the dipole moment in the respective 
component direction. 

Battan [47] ".efines the dipo<e moments in terms of particle shape 
as follows: 


= 41- .^tg. g'}E^ 

m 

= 4ir £^{g, g'lE^ 


(6.1-19) 


where g is i’ complex number describing the effects of particle 
geometry on the dipole moment associated with the particle's axis of 
revolution and g' is associated with the orthogonal axis. The defini- 
tions of g and g‘ are given below: 


4ir (e^ - 1)P 

g* = ~ 

4ir + (e^ - 1)(2 it - P/2) 
where for a prolate spheroid, 


( 6 . 1 - 20 ) 


( 6 . 1 - 21 ) 






Ice Particle 



Figure 6.1-2. 




I 


iinhi 



Geometry used in the derivation 
of the single ice particle 
scattering coefficients. 
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’’pro " ‘ ( 6 .!- 22 ) 

for an oblate spheroid, 

W ■ ^ (6.1-23) 

6 e 

and for a spherical particle, 

P ■ y . (6.1-24) 

The various symbols appearing above are defined below: 

complex refractive index of air (s^q 3 1 } 

: complex refractive index of ice [47] 

(cp ■ 1.73 - j 0.0024 $ 0® C}* 

Vp : volume of the ice particle 

e : eccentricity fe » 1 * b^/a^} . 

the eccentricity is close to one, i.e., a long narrow prolate 
spheroid or a flat platelike oblate spheroid Eq. (6.1-22) and Eq. 
(6.1-23) reduce to 

PpRO ■ ® 

(6.1 -2S) 

’’OBL * • 

Using Eq. (6.1-24) and Eq. (6.1-25) In Eq. (6.1-20) and Eq. (6.1-21) 

* Haworth, McEwan, and Watson [30] evidently have a misprint in their 
publication. 






■ 'I'nnaMw 
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and th« results In Eq, (6J«19) yield for the various particle shapes; 


Prolate 




h ’ I 


d '■ h 


7 (4 - 1) M 

2 -4 4it e 

(4^y) ^ ^ 


^ * f ^ ,-2 


Oblate 


I«h • - n 4 m E? 




^0 "h 


"v * ‘o 


Sphtrlcil 


"h.y • 


(«R + i) 


4n E 


0 S.V 


where 


(6.1-26) 


(6.1-27) 


(6.1-28) 


¥ 1s the equi volume trie particle radius In 

Now consider the radiation properties of an electric dipole. From 
Stratton [48] the radiated electric field Is 

«RAO ki ^ . '^V (6,1-29) 

E(j) • [l X (m K l)] — j — 

*0 
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nRAD . '^0 „ e ° 

ITT % “T“ 


.RAD 


° m. ® 


(6.1-30/ 


4ir "*v 2 

0 


Comparing Eq. (6,1-30} to the definition of the scattered field as 
defined by the scattering model, the single particle scattering 
coefficients can be determined as seen below: 


^h,v * ^h.v ^h.v 2 


-jk 2 

® = (z) 

‘^h.v 


-jk^z 


nil, 


4u Eq h,v 2 


So. 


h.v 


4tt a. 


"h.v 


(6.1-31) 


Using Eqs. (6.1-26). (6.1-27). and (6.1-28) in (6.1-31) yield the 
single particle scattering coefficients for each particle shape and 
are given below: 


Prolate 

. i±L-p> (6.812211 - j 0.026842} 

(6.1-32) 

/RO ^ L4tt P} (3^268519 - j 0.006179} 

^ 3r 

Oblate 

Ml _ {4tt a^} 

I’ 


{6.812211 - j 0.026842} 
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Spherical 


fOBL « {2.150047 - j 0.002674) 

^ 3X 


fSPH ^ a ) { 3 , 95417 . j 0.009044} . 
3X^ 


Equations (6.1-32), (6.1-33), and (6.1-34) define the single ice 
particle scattering coefficients. 8y summing over all the particle 
radii within an ice particle concentration, a general scattering 
coefficient can be obtained. Define 


''h.v ° _ 

a = a. 


n(a) 


(6.1-35) 


where n(¥) is the number of particles with radius a. Since fu is 

n jv 

directly proportional to the volume of the particular particle size 


class and inversely proportional to 


F s A 
h,v h.v 


(6.1-36) 


where V is the Uital volume of Ice in a cubic meter of air and Au .. 

n,v 

is a complex constant Independent of frequency given in Table 3.6-1. 
This result along with the general scattering coefficient formulation, 
Eq. (3.1-1), and the elevation angle formulation presented in 3.6 yield 
Eq. (3.6-12). 


6.? The Rain Propagation Prediction Program 
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6,2.1 General Description 

The Sain FnpapagaHcn prediaHon pi^grm uses the scattering model 
and the synthetic storm algorithm to predict attenuation, isolation, 
and phase as a function of rain rate for a given set of syst^ and 
rain parameters. A block diagram of the RPP program is found in Fig. 

6.2- 1. The single input format of Fig. 6.2-2 facilitates changing the 
system and rain parameters. Although the lack of data has caused 
certain restrictions to have been placed on the particle distribu- 
tions along the rain extent, the scattering model and thus the RPP 
program can model an arbitrary inhomogeneous rain given the appropriate 
input parameters. The input parameters will be discussed in more 
detail in Sec. 6.2-2. 

After reading the input parameters and calculating the clear 
weather isolation and phase, the RPP program prints an information 
header at the beginning of the output section as seen in Fig. 6.2-3. 

The rain conditions for each cell are then computed and the scattering 
matricies as defined by Eq. (3.3-1) are determined. After computing 
the wave exiting the rain medium and the wave-antenna interaction, the 
RPP program outputs rain rate, attenuation. Isolation and phase data 
as seen In Fig. 6.2-4. Since curve fits to theoretical data are often 
useful, the RPP program outputs curve fit data for various combina- 
tions of rain rate, attenuation, and isolation data as seen in Fig. 

6.2- 5. 


OF rOOTl 



QRIGINAL 

OS fiofift auwJiafi 


Figure 6.2-1. Block diagram of 

the rain propagation 
prediction program (RPP). 
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IKO 

33.0 

1000010000 1000 
10 10 10 

3 60 3 

-43.35 0.00 

-44.01 0.00 44,11 90.00 

0.0 0 . 0 - . 66 - . 66 - . 66 - . 66 - . 66 - . 66 - . 66- . 66 
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 
0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 
12. 012.012. 0,:.012.012. 012. 012. 012.012.0 
0000000000 
0 



Figure 6.2-2. Sample input data to the rain 
propagation prediction program 
(RPP). 
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STOP 
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10 


10 
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J 


60 
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OlNKEO COEFFICIENTS 









T 
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4 
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10 

RAIN RATE 

0.0 
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12.0 
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O 

0 

O 

0 
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0 

0 

0 
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• 44v« CLEAR U£Atrit:H IrlKiRHAllUH *«*•« 
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THE fULAHllAlllIH PAHAULTEK EPSILUI PIIR 7riE Of ANTENNA PIILAK 1 2 AT ItIN STATE iS 

THE PHLARI2ATIUN PArtAHLIER TAU EuK Tub CU ANIENNA PUL ARAZA1 I(iH STATE IS 

THE POLARIZATION PARAKtfbH EPSILON FOR THE UlUSS ANTENNA POLARIZATION S^Alb IS 

THE POLARIZATION PARAhETtH TAU EUR Tilt CROSS ANTENNA SlATt IS 

THE CLtAK HEAIMLK ISULATIUN RtFtRbNCllI TO IHb ANTENNA PURlS IS 

THE CLfcAH IJEAlllbH PHASE UIEFLKtHLE 
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0.0 OEGRbtS 
-44.01 UEGREES 
0. 0 OEbREE S 

44.11 lltGREtS 
VO. 00 OEGREEs 
37. » UU 
-C.OI OEGPtbS 


Figure 6.2-3. Information header section for the output of the 
rain propagation prediction program corresponding 
to the Input data of Fig. 6.2-2. 
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aTTI'N fOul 

UOL <0BI 


PHASk lUfcGl 
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1.93 

31.03 


299.96 

b.OO 

2.90 

29.51 


290.61 

9.00 
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285.01 

12.00 
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26.63 


283.15 
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24.60 


279.79 

21.00 

6.21 

23.94 


270.72 

24.00 
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23.24 
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27.00 

7.26 

22.66 


276*90 

30.00 
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22.12 


276.30 
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275.65 
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39.00 
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20.66 

. 
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42.00 

9.61 

20*20 
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43.00 

10.07 

19.76 


273.60 

tO.OO 

iu.$l 

19.37 
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31.00 

10.94 

19.01 
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34.00 

11.34 

18 .68 


272. 5U 

57.00 . 

11.77 

18.32 
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bO.OO 

12.20 

17.90 
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Figure 6.2^4. Kain propagation prediction program output, attenuation. 
Isolation and phase as a function of rain rate 
corresponding to the input data of Fig. 6.2'2. 
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Figure 6.2-5. Rain propagation prediction program output, 
curve fits to the data In Fig. 6.2-4. 
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6,2,2 RPP Users Guide 

The RPP program was designed to provide an efficient tool to 
predict the effects of precipitation on satellite comnunl cations Siystem 
performance. The RPP program Is written In the FORTRAN programming 
language, and the program listing is found in Sec, 5.2,3. This 
section explains the steps In the operational procedure of the RPP 
program. Each step corresponds to a line of input data In Fig. 6.2-2, 

Step 1 

Input the operat1r.g frequency in GHz in an F5.1 format.* 

Step 2 

Input the system elevation angle in degrees in an F5.1 
format. 

Step 3 

Input the effective physical rain extent (L) do loop 
parameters in meters, LSTET^ LSTOP^ and ljNC* The input 
format is 15. The do loop parameters allow finer control 
over the rain extent parameter. To implement the synthetic 
storm algorithm consult Fig. 3.6-5, 


* The scattering coefficients f^ and arc limited to 11,0, 14.0, 
20.0, and 30.0 GHz. To increase the frequency selection, curve 
regressions for published fy and fh coefficients at other frequen- 
ies are easily obtained and the RPP program is easily modified. 
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Step 4 

Input the mean canting angle (?^) d locfg parameters AUSTET^ 
ANSTOF, and MINC in degrees {|?J < 45®}. The format is 15, 

Step 5 

Input the rain rate do loop parameters RRSTST^ SRSTOP^ and 
BRINC in mm/hr {RR < 250 mm/hr). The format is 15. 

Step 6 

Input the clear weather wave polarization state parameters 
and in degrees. The format is F7.2, 

Step 7 

Input the CO and cross antenna polarization state parameters 
^CO- ^co* "cross "cross degrees. The format Is F7.2. 

Step 8 

Input the X.. coefficients of the synthetic storm algorithm 
for each rain cell. To model a uniform rain path * 0. 

Note that aetl one is aZoseet to the reaeive antenna* The 
format is F4.1. 

Step 9 

If a distribution of mean canting angles is desired along 
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the propagation path, input che multiplicative factors for 

each T>avn cell {? = If a uniform distribution 

^i ■ ^1 

along the rain path is desired = 1.0. The format is 

F4.1. 


Step 10 

Input the frastion of oblate raindrops in each rain 
cell. The format is F4.1. 

Step 11 

Input the standard deviation of the canting angle distribu- 
tion of a- of each i>ain cell in degrees. The format is F4-1. 

Step 12 

Input a "0“ to indicate a rain cell or input a "1” to 
indicate an ice cell for each cell along the propagation 
path. For an ice cell, steps 9, 10 and 11 are overridden 
and the assumptions of Sec. 3.6.5 apply. The format is 13. 

Step 13 


Input a "0" for a downlink. Input a "1" for an uplink. The 
format is 15. 

By controlling steps 1 through 13, a wide range of system and 
rain parameters can be modeled and their effects on the performance 
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of an arbitrary satellite communications link can be determined. 


T 

) 


c 

C 6.2.3 RPP Program Listing 

C 

c* 

C* ;:iAIN PROPAGATION PREDICTION PROGRAM IRPPI 

C* 

C* THIS PROGRAM OUTPUTS THEORETICAL RAIN PROPAGATION 
C* DATA {ATTENUATION. ISOLATION. AND PHASE) FOR A 

VARIETY OF FREQUENCIES, ELEVATION ANGLES AND LOCATIONS 
C* AS A FUNCTION OF GROUND RAIN RATE. 

C« THE INPUT FORMAT ALLOWS THE RAIN PATH TO BE SUBDIVIDED 

C* INTO TEN SEGMENTS WHICH PROVIDE ARBITRARY VARIATION 

C* OF THE PROPAGATION MEDIUM PARAMETERS- 

C* 

C* PROPAGATICT,' MODEL: SCATTERING MODEL 

C* 

C* PROGRAMMER: RUSSELL R« PERSINOER 

C* 

C* CONSULTING ADVISORS: OR* WARREN L* STUTZHAN 
C* AND ROBERT E* CASTLE 

C* 


RPPOOOlO 

RPP00020 

RPP00030 

RPPOOOAO 

RPP00050 

RPP00060 

RPP00070 

RPP00080 

RPP00090 

RPPOOlOO 

RPPOOIIO 

RPP00120 

RPP00130 

RPPOOIAO 

RPP00I50 

RPP00160 

RPP00170 

RPP00I80 

RPP00I90 

RPP00200 

RPP00210 

RPP00220 


C* THIS PROGRAM WAS DEVELOPED IN PARTIAL 
C* FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE 
C* OF MASTER OF SCIENCE IN ELECTRICAL ENGINEERING; 

C* HOWEVER, ITS DEVELOPMENT WAS SUPPORTED BY NASA 
C* WITH important ADDITIVE SUPPORT FROM THE DEFENSE 
C* COMMUNICATIONS AGENCY AND THE US ARMY RESEARCH OFFICE, 
C* 

Q4^^iiti:tt**)^***tti4i***********4^*^*4‘*^^**************>^************* 

C 

C 


RPP00230 

RPP00240 

RPP00250 

RPP00260 

RPP00270 

RPP00280 

RPP00290 

RPP00300 

RPP00310 

RPF00320 



r> o 


C 


DIMENSION ICEf lOl.SIGHAIlO) 


RPP00330 

DIMENSION OIN(iO) 


RPP00340 

DIMENSION AllOJ fB( 10),CI 10! 


RPP00350 

DIMENSION THETA(15)»RRATE(15l,RRATEU300)«PEn5} 

RPP00360 

DIMENSION PHASEl(300l«IS0Llf300}»ATTENlI300) 


RPP00370 

DIMENSION PXXSlSIfPXYl 151 tPTXf 151 »PYY( 151 


RPP00380 

INTEGER LSTRTtLSTOP«LlNC 


RPP00390 

INTEGER ANSTRTfANSTOPfANlNC 


RPP00400 

INTEGER RRSTRT»RRSTOP«RRINC 


RPP00410 

INTEGER BINtOlREC 


RPP00420 

REAL I SOL 


RPP00430 

REAL ISQLltLOKANG 

BPP00440 

COMPLEX MINUSJ 

RPP00450 

C0MPLEX4C16 CONSNT 

RPP00460 

C0MPLEX«16 EXO«EYOtEXCvEYC»EXX«EYX 

RPP00470 

COHPLEX^iS FV*FH«FSPH«DIFF 


RPP004B0 

C0MPLEX»16 PI»P2«P3tP4 


RPP00490 

COMPLEX* 16 EXNfEVN 


RPP00500 

C0MPL£X*16 PXX,PXV,PYX.PVY 


RPP00510 

COMPLEX* 16 <^XX,FXY*FVY 


RPP00520 

C0MPLEX*16 SXXtSXYtSYXtSVV 


RPP00530 

COMPLEX* 16 LAMOAl, LAHDA2 ^ALPHAOe ALPHAl 


RPP00540 

RPP00550 

COHHON/BLOCl/MINUSJ.CONV 


RPP00560 

RPP00570 

RPP00580 




C* 

C* fi:*sr THE PROGRAM READS ALL NECESSARY VARJABLE VALUES INTO 

C« HEHORY AND WRITES THIS INFORMATION OUT AS A HEADER AT THE 

C« BEGINNING OF THE OUTPUT SECTION. 

C* ALSO INCLUDED IN THIS SECTION ARE THE CLEAR WEATHER CONDITIONS 


RPP00600 

RPP00610 

RPP00<^20 

RPP00630 

RPP00690 



c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

€♦ 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


AND AS A RESULT, UNIT VECTORS DESCRIBING THE SYSTEM POLARIZATIONS RPP00650 
ARE COMPUTED, RPP00660 


FREQ IS FREQUENCY IN GHZ ( 11,0, 14.0, 20.0, 30.0 ONLY I 
LOKANG IS ELEVATION ANGLE IN DEGREES CO, ..901 
LSTRT,LSTOP, AND LlNC ARE THE RAIN EXTENT PARAMETERS 
IN METERS 

ANSTRT, ANSTQP, AND ANINC ARE CANTING ANGLE PARAMETERS 


RPP00670 

RPP0C680 

RPP00690 

RPP00700 

RPP00710 

RPP00720 


IN DEGREES 

RRSTRT, RRSTOP, AND RRINC ARE RAIN RATE PARAMETERS IN MM/HR 

ALL EPS'S AND TAU« S ARE IN DEGREES 

BIN ONE IS LOCATED AT THE GROUND ANTENNA 

All,. ..101, B(l,...101, AND Cfl,...10) ARE MULTIPLICATIVE FACTORS 
ASSOCIATED WITH THE RAIN EXTENT BINS 
CONTROLS THE ITH BIN RAIN RATE 


All I 

THE VALUES OF A(U ARE SET IN ACCORDANCE WITH THE 


rtrKuu f w 
RPP00750 
RPP0076C 
RPP00770 
RPP00760 
RPP00790 
STORM ALGORITHNRPP00800 


BID CONTROLS THE ITH BIN CANTING ANGLE 


RPPOOBIO 


cm CONTROLS THE ITH BIN PERCENT 06LATENESS 

THE ABSOLUTE VALUE OF THETAII) IS ALWAYS LESS THAN 45 DEGREES 
RRATEdl IS ALWAYS LESS THAN 250 MM/HR 
SIGMA! !) IS THE STANDARD DEVIATION OF THE CANTING ANGLE 
DISTRIBUTION FOR THE ITH BIN. lEQUAL TO ZERO OR 10<S1GMA<50I 
IF ICE(I) IS EQUAL TO 1, THE ITH BIN IS FILLED WITH ICE PARTICLES 
TF ICEin IS EQUAL TO 0, THE ITH BIN IS FILLED WITH RAIN 
OIREC EQUALS ONE F'^R UPLINK 
DIREC EQUALS ZERO FOR DOWNLINK 


RPP00820 

RPP00830 

RPP00840 

RPP00850 

RPP0C86O 

RPP00870 

RPP00660 

RPP00890 

RPP009D0 


RPP00910 




C 


READ! 5, 100) 
R£A0I5,100| 
REA0I5,10L) 


FREQ 

LOKANG 

LSTRT,LST0P,L1NC 


RPP00930 

RPP00940 

RPP00950 

RPP00960 


REA0I5 
READ! 5 
READfS 
READ! 5 
READ (5 
READ I 5 
REA0(5 
READ 1 5 
REAO(? 
REA0I5 

100 FORMAT 

101 FORMAT 

102 FORMAT 

103 FORMAT 

104 FORMAT 

105 FORMAT 


1011 ANSTRTtAr^STOPf AMINC 

101) RRSTRT,RRST0P«RR1NC 

102) EPSWtTAUW 

103) EPSCfTAUCtEPSXtTAUX 

104) (A<n»l»lvl0) 

104) |8(1)»I==1«10) 

104) (C( Dtlsl.lO) 

104) CS1GMUI)«1«1»10) 

105) lICEfI ),1=1»10) 

1011 OIREC 

F4* 1) 

ISf IS|I5) 

F7.2,F7.2) 

4(F7,2) ) 

10(F4.in 

1013) 


WRITE (6,200) 

WRJTE(6,221) 

IF(OIKEC.EQ.I) WRITE(6,203U FREQ 
IF(01REC«EQ.O) URITE(6t2030) FREQ 
URlTE(6f201) LOKANG 
WRITE(6,221) 

WRITE(6,202) 

WR1TE(6,204) 

WRITE (6,205) 

WR1TE(6,206) LSTRT ,LSTOP ,LlNC 
WRITE (6,207) ANSTRT, ANSTQP , AMINC 
WR1TE(6,208) HRSTRT,RRSTOP,RRINC 
WRITE(6,20<9) 

WRJTE(6,210) 



RPP00970 

RPP00980 

RPP00990 

RPPOIOOO 

RPPOlOlO 

RPP01020 

RPP01030 

RPP01040 

RPP01050 

RPP01060 

RPP01070 

Rppoioeo 

RPP01090 

RPPOllOO 

RPPOillO 

RPP01120 

RPP01130 

RPP01140 

RFP01150 

RPP01160 

RPP01170 

RPP01180 

RPP01190 

RPP01200 

RPP01210 

RPP01220 

RPP01230 

RPP01240 

RPP01250 

RPP01260 

RPP01270 

RPP012&0 



WR1TEC6«2UI IA(ll,I»l9l01 
WR1TE(6»212} (Bflltl^lrlOl 
WRlTEi6»213) (C(ll vl^lflOl 
WRITEI6.2131I f SlGHAClItl^lflOl 
WfllTEC6t2132) (ICEIIl, 
WR1TE(6,214I 
WRlTE(6v215) EPSW 
MR1TEI692I6I TAUU 
URITEI69217) EPSC 
WRITE (69218) TAUC 
WR1TE(6»219) EPSX 
WR1TE(69220) TAUX 


RPP0129D 

RPP0I300 

RPP01310 

RPP01320 

RPP01330 

RPP01340 

RPP01350 

RPP01360 

RPPQ1370 

RPP01380 

RPP01390 

RPP01400 

RPP01410 


200 FORMAT aHlt3X, *»♦♦♦♦ RAIN PROPAGATION PREDICTION PROGRAM CRPP) ♦♦♦RPP01420 

) RPP01430 

201 FORMAT (/»3X« "THE ELEVATION ANGLEf BETA, IS • .FA. 1«2X« •DEGREES* ) RPP01440 

202 FORMAT (// 9 3 X 9 PROGRAM PARAMETERS ♦♦♦♦♦•) RPP01450 

2031 FORMAT I //93Xf* THE UPLINK PROPAGATION FREQUENCY 15«9 2X9F4*l92X«*6HZRPP01460 


l*l 


INCREMENT* I 


2030 FORMAT I //9 3 X 9 * THE DOWNLINK PROPAGATION FREQUENCY IS *92X»F4*l92X9 
1*GHZ* ) 

204 F0RNATI//t3X«*OO LOOP PARAMETERS*) 

205 FORHATI/93X,* TYPE START STOP 

206 F0RMATI/93X9*!) RAIN EXTNT * 9l695Xtl694X9l6) 

207 F0RMATI/93X9*2) CANT ANGLE * 9l695Xtl694Xfl6) 

208 FORMAT (/ 9 3X««3) RAIN RATE * 91695 X 91694 X 916 ) 

209 FORMATI/Zt 3 X 9 * BINNED COEFFICIENTS*) 

210 F0RMAT(//9l5X**l*,5X,*2*f5Xt*3*,5X,*4*95X,*5*#5X9*6*95X,*7*95X9*8*RPP01560 

195 X 9 * 9 * 95 X 9 * 10 *) RPP01570 

211 FORMAT (/• 3 X»*R AIN RATE* , 2 X 9 101 F4«l92X) ) RPP01580 

212 FORMATI/.SXv'CANT ANGL* * 2 X 9 lOf F4. I 92 X) ) RPP01590 

213 FORMAT!/ 93Xt*FRAC OBLT* 92X«10(F4.1«2X)I RPP01600 


RPP01470 

RPP01480 

RPP01490 

RPP01500 

RPP01510 

RPP01520 

RPP01530 

RPP01540 

RPP01550 


2131 F3RI4AT(/»3X»«S1GHA 

2132 FORMAI (/«3XtMCE 
21A FORMAT(//»3X, 


■ ,2X,10(F4.1»2X)1 
•«2X»10(14.2XI ) 

CLEAR UEATHER INFORHATION 




215 FORMAT C//»3X, •THE POURIZATION PARAMETER EPSILON FOR THE INCIDENT 


IHAVE IS* »21X,F7.2,2X,>DEGREE5* I 

216 FORHATI/p3X«*TH£ POLARIZATION PARAMETER 
1 IS*. 24X»F7.2»2Xt*OEGREES« ) 

217 FORMAT {/«3X,* THE POLARIZATION PARAMETER EPSILON 
1 POLARIZATION STATE IS • » 5X.F 7,2 .2X, 'DEGREE S* I 

218 FORMATI/.3X,*THE POLARIZATION PARAMETER TAU FOR 
lARAZATlON STATE IS* ,9X,F7.2.2X,*DEGREES* 1 

219 FURMATt/.3X»*THE POLARIZATION PARAMETER EPSILON FOR 
INNA POLARIZATION STATE 1 S* . 2X, F7,2 ,2X, *DEGRtES * I 


RPP01610 
RPP01620 
RPP01630 
RPP01640 
RPP01650 

TAU FOR THE INCIDENT HAVE RPP01660 

RPP01670 

FOR THE CO ANTENNARPP0I680 

RPP01690 

THE CO ANTENNA P0LRPP01700 

RPP01710 
THE CROSS ANTERPP01720 
RPP01730 


220 FORMAT I /. 3X. * THE POLARIZATION PARAMETER TAU FOR THE CROSS ANTENNA RPP01740 

I STATE IS* ,l8X,FT,2t2X, •DEGREES*I RPP01750 

221 FORMAT ( //. 3X . * ^lf^*^*********‘*:^^***************^**^********************RPPQ 1 760 

I RPP01770 

222 FORMAT! 1H1.2X.* THE RAIN EXTENT IS *.I5#* METERS* » RPP01780 

223 F0RHATI//*3X,*THE MEAN ANGLE BETWEEN THE Y-AXIS AND THE DROP MIN0RRPP01790 

I AXIS ITHETAJ IS *,I5.* DEGRESS*! RPP01800 

SCATTERING MODEL ♦♦♦»♦♦♦♦♦*♦♦♦♦♦*, RPPOISIO 

RFPO102O 
fDBt • .8XRPP01B30 
RPP01840 
RPP01850 
RPP01B60 


224 FORMAT I////, 25X, •♦*♦**♦★♦**♦♦♦ 

I/I 

225 FORMAT!//, 3X.*RAIN RATE !MM/KR I • , 5X, » ATTEN I OB 1 • ,9X. • I SOL 
1. 'PHASE !OEGI*,//l 

226 FORMAT!/. I0X.F6, 2. IIX. F6.2. 12X.F6, 2. tlX.F7.2) 

227 FORMAT!//////* CURVE FITS TO ABOVE DATA* I 


228 FORMAT!///, 21X, 'CURVE 
I* ,9X, *R**2* 1 

229 FORMAT!//, • ATTEN VS 
18,4,4X11 

230 FORMAT!//,* ISOL VS 
18.4,4X1) 


TYPE* ,6X,*ASS0CiATE0 EQUATION* ,7X,*A* ,11X, *BRPP01870 

RPP01860 

RR* ,8X, 'POWER*, 12X, *ATTEN - A^BR**B* ,4X,3!FRPPO109O 

RPP01900 

RR'.lOX, 'LOG*. 12X,*IS0L = A+B*LNIPR) * ,3X, 3! FRPP01910 

RPP01920 


I 

CT* 

ro 

■ 


231 FORHATI//, • 
lXt3(F6.4«4X)| 


ISOL VS ATTEN*t?X«*tOGSlOX, •ISOi. - 


A4^B«LNCATTENI S3RPP01930 

RPP01940 

C RPP01950 

C • RPP01960 

RPP01970 
RPP01960 
RPP01990 
RPP02000 
RPP02010 
RPP02020 
RPP02030 
RPP02040 
RPP02030 
RPP02060 
RPP02070 
RPP02080 
RPP02090 

C RPP02100 


PI=3. 14159265 
CONV-P 1/180. 

LOKANG=LOK ANG«CONV 
Mlf4U5J-tO« «-l.) 

CALL C0HPNTfEPSM,TAUH«EX0«EY01 
CALL C0MPNT|EPSC,TAUC,EXC*EYC1 
CALL CaNPNT(EP$X>TAUX«EXX«EYX) 

EYC'OCONJGfEYC} 

EYX-0C0NJ6fEYXI 

CALL CLEW£ACEXO»EYO»EXC«EYCfEXX»£YXt 


C* 

C* 

C» 

C* 

C* 


THE REST OF THE PROGRAM USES THE SCATTERING MODEL 
TO COMPUTE ATTENUATfONt ISOLATION* AND PHASE AS A FUNCTION 
OF RAIN RATE AND OUTPUTS THEM IN AN ORGANIZED FORM. 


RPP02120 

RPP02130 

RPP02140 

RPP02150 

RPP02160 


C RPP02180 

RPP02190 
RPP02200 
RPP02210 
RPP02220 

+ AN INC ’ RPP02230 

4 - AN INC RPP02240 


IFIANSTRT.LT.OI 

IFIANSTRT.LT.O) 

IFIANSTRT.LT.OI 


ANSTRTsANSTRT 

ANSTOP®ANSTOP 


FLAG«l.O 

ANSTOPsANSTOP4^360 

ANSTRT*ANSTRT+360 

AN INC 
AN INC 



Oi 

CO 


non non 


C 

00 I ll^LSTftT.LSTOPtLlNC 
00 I JJ-ANSTRT,ANSTOP»ANfNC 
C 

RArNL-II 

jjj=Jj-ANINC 

IFfFLAG*£Q.1.0| JJJJ=JJJ-360 
C 

WRITE(6,222I 11 
WRITE (6,223) JJJJ 
WR1TE(6,224) 

WR1TE(6,225) 

C 

Ill-O 

THETA (U®JJJ»CONV 

IF(THETAm.GT.Pl) THETA (1 ) =THETAI l)-2.*PI 
C 

00 2 KK»RRSTRT,RRSTOP,RRfNC 
SET UP THE BINNED PARAMETERS 
RRATE(I)»KK 


C 

C 


III IS THE NUMBER OF RAIN RATE ENTRIES... LESS THAN 300 

!1I>1I1 * I 
RRATEUlll )==KK 

00 3 !al,IO 

K=l + 1 


RPP02250 
RPP02260 
RPP02270 
RPP02280 
RPP02290 
RPP02300 
RPP02310 
RPP02320 
RPP02330 
RPP02340 
RPP02350 
RPP02360 
RPP02370 
RPP02380 
RPP02390 
BPP02400 
RPP02410 
RPP02420 
RPP02430 
RPP02440 
RPP02450 
RPP02460 
RPP02470 
RPPO240O 
RPP 02490 
RPP02500 
RPP02510 
RPP02520 
RPP02530 
RPP02540 
RPP02550 
PPP02560 



non n o n o o non 


IFlK.EQ.llt GO TO 9 RPP02570 

RRAT£(IC}-4ftRATEm/lO.O)4c*A<K)«RRATE(l) RPP02580 


THETAIK)«BIKI*THETA<11 RPP02590 

9 PenJ-Cdl RPP02600 

C RPP026L0 

IFIRRATEIKI.LT.1.00) RRATEIK)-1,0 RPP02620 

IFtRRATECKKGT. 250.0) RRATE C K) =250.0 RPP02630 

IF * THETA CKI.GT. 0.7854) THETA IK »«0. 7854 RPP02640 


I F C THETA I K) . LT .-0. 7854 » THETA I K) =-0. 7854 
3 CONTINUE 

SCATTERING MODEL 

OELTAL=l.O 

C0NSNT=MINUSJ*DELTAL*0.3/FREQ 
NN=RA1NL/IOELTAL’*'10.0) 

COMPUTE THE SYSTEM MATRICES FOR EACH DISCRETE BIN 
USING THE CAYLEY-HAMILTON THEOREM, 

11=0 

BINm = L 
DO 4 I=l»10 

CHECK FOR UNIQUE BINS, IF PREVIOUSLY COMPUTED DO NOT RECOMPUTE 

IFU.EQ.l) GO TO 13 
C 

DO 8 K1=1,I1 
C 

K2=BIN{KL) 


RPP02650 

RPP02660 

RPP02670 

RPP02680 

RPP02690 

RPP02700 

RPP02710 

RPP02720 

RPP02730 

RPP02740 

RPP02750 

RPP02760 

RPP02770 

RPP02780 

RPP02790 

RPP02B0O 

RPP02810 

RPP02820 

RPP02830 

RPP02840 

RPP02850 

RPP02860 

RPP02870 

RPP02880 


CTl 

U1 


RPP02890 

IFU THETAC It.EQ.THETAC K2n.AND.fPEm.EQ.PE(K2)l.AND.(RRATEf n RPP02900 

1*EQ.RRATE(K21I .AND. nCES n.EQ. ICE (K2) I *ANO. (SIGMA (I).EQ. SIGMA fK2)IRPP029l0 
1] GO TO 7 RPP02920 

RPP02930 


6 CONTINUE 

13 CONTINUE 

AM0OE=0.5 i-O.ASOALOGHHRRATEd n 
IF(1CE( n.NE.l) GO TO 14 

CALL lCECOF(FREQrANODE»FXXiFYY)FXYfLOKANG»RRATE»I »THETA) 

60 TO 15 

14 CALL CAPF(PEiAHOOEtFXXiFYYfFXY«FREQ»LOKANCfRRATEf ItTHETA.SIGHAJ 

15 CONTINUE 

SX X= I • ♦r.ONSNT*FXX 
SYY=l.*CONSNT*FYY 
SXV=CONSNT*FXY 
SYX*SXY 

LAMDAl=ISXX«YYl/2.+C0SQRTJ I (SXX+SYYI**2I/4.4-SXY+*'2~SXX*SYY> 
LAMDA2*(SXX+SVYJ/2 .-COSQRTH (SXX+SYY )*«=2 >/4.4^SXY**2-SXX*SYY) 

ALPHA 0-( LAM0A2* e LAMDAl ♦*NN ) -LAMDAl LAM0A2**NN »» / ( LAM0A2-LAMDA1 1 
ALPHA I^ILAHOA^'^^'NN - L AHDA1**NNS / 1 LAMDA2-LAM0AU 


PXXU )=ALPHAO ALPHAl^SXX 

PXV(n = ALPHA' SXY jQ q 

PYX(n-pxY(n 

PVYin*ALPHAO+ALPHAl»SYY SC 


RPP02940 

RPP02950 

RPP02960 

RPP02970 

RPP029eo 

RPP02990 

RPP03000 

RPP030I0 

RPP03020 

RPP03030 

RPP03040 

RPP03050 

RPP03060 

RPP03070 

RPP03080 

RPP03090 

RPP03100 

RPP02U0 

RPP03120 

RPP03130 

RPP03140 

RPP03150 

RPP03160 

RPP03I70 

RPP031B0 

RP>^03190 

RPP03200 



11 = 11+1 
B1NC11)»I 
GQ TO 4 
7 CONTINUE 
C 

PXXm=PXXIK21 
PXVm=PXYIK2J 
PYX( I )=PYXIK2I 
PYYU J*PYY<K2I 
C 

4 CONTINUE 
C 

DO 5 I=lt9 
C 

IFCOlREG.EQ.ll GQ TO 11 
C 

J=ll-I 
K=J-1 
GO TO 12 

11 CONTINUE 
K=I+l 
J=I 

12 CONTINUE 
C 

Pl>PXXfK34<PXXfJl + PYX(K)«PXY( J) 
P2=PKX(K)*PVXI J) + PYXIKJ*PYVU) 
P3-PXYfK|4tPXXlJl + PYY(KI*PYXIJ1 
P4=PXYIKI*PYXI J) + PVyiK>*PYYIJI 
C 

PXXIKI^PI 

PYX(KI=P2 



RPP03210 

RPP03220 

RPP03230 

RPP03240 

RPP03250 

RPP03260 

RPP03270 

RPP03280 

RPP03290 

RPP03300 

RPP03310 

RPP03320 

RPP03330 

RPP03340 

RPP03350 

RPP03360 

RPP03370 

RPP033B0 

RPP03390 

RPP03400 

RPP03410 

RPP03420 

RPP03430 

RPPD3^,40 

RPP03450 

RPP03460 

RPP03470 

RPP034BO 

RPP03490 

RPP03500 

RPP03510 

RPP03520 



r»nn non n non non 


C 


PXYCK)=P3 RPP03530 

PYY(KJ=P4 RPP03540 

RPP03550 

5 CONTINUE RPP03560 

RPP03570 

COMPUTE THE RESULTJiNG X AND Y COMPONENTS EXITING THE RAIN MEDIUM RPP03580 

RPP03590 


EXN=Pl*eXO + P2»EY0 RPP03600 

EYN=P3*£X0 + P4*EY0 RPP03610 

RPP03620 

COMPUTE THE ANTENNA OUTPUTS DUE TO RAIN AND POLARIZATION MISMATCH RPP03630 

RPP03640 

CALL OUTANT(EXNfeYN#EXC,£YC*EXX,eYX,III ,IS0L1,ATTEN1,PHASEU RPP03650 

RPP03660 

2 CONTINUE RPP03670 

PPP03680 

NOW OUTPUT THE DATA IN AN ORGANIZED FORM RPP03690 


DO 6 1=1,111 

WR1TEC6,226) RRATEl U) ,ATTEN1 f II ,ISOLU I) , PHASEll I) 
6 CONTINUE 

PERFORM LINEAR REGRESSIONS ON OUTPUT DATA 

WRIT£<6,227) 

WR1T£(6,228) 

CALL CURVE (RRATEL, ATTENlfl !I«AA,BB,R, 21 
WRITE<6,229| AA,BB,R 

CALL CURVE (RRATEl, ISOL L , II I , AA tBB,R,3I 
WRITE(6,230I AA,8B,R 

CALL CURVE I ATT ENl, IS0L1,11 1 ,AA,66«R,3) 

WRITE(6,231I AA,8B,R 


RPF03700 

RPP03710 

RPP03720 

RPP03730 

RPP03740 

RPP03750 

RPP03760 

RPP03770 

RPP03780 

RPP03790 

RPP03800 

RPP03810 

RPP03820 

RPP03830 

RPP03840 


I 


I 


168 



o n o o 


C 

I CONTINUE 
C 

MR1TE16>IOOO) 

1000 FORMAT (IHll 
STOP 
END 

SUBROUTINE COMPNT(EPS,TAU,EX,EY) 

THIS SUBROUTINE RETURNS THE X AND V COMPONENTS GIVEN AN EPSILON 
AND TAU UN DEGREES) DESCRIBING AN ARBITRARY POLARIZATION STATE 

C0MPLEX*16 EX,EY 
COMPLEX HINUSJ 
COHMON/BLDCl /MINUS JtCONV 
C 

EP5R=EPS*C0NV 

TAUR=TAU*CONV 

C 

IFIABSlEPSR»-EQ.(45.*C0NVn GOTO 1 
IFIEPSR.EQ.O.) GO TO 2 
IFCTAUR-EQ.O,) GO TO 3 
!F(TAUR.EQ.{90.*C0NV)J GO TO 4 
C 

Tl=TAN(2.*EPSR) 

T2=S1NI2.*TAUR » 

DELTR^=ATAN2ITl*T2) 

GAMR=0.5*ARC0S(C0S(2*EPSRJ*C0St2*TAURJ ) 

GO TO 100 
1 DELTR=2.«=EPSR 
GAMR=45**C0NV 
GO TO 100 


RPP03850 

RPP03860 

RPP03870 

RPP03880 

RPP03890 

RPP03900 

RPP03910 

RPP03920 

RPP03930 

RPP03940 

RPP03950 

RPP03960 

RPP03970 

RPP03980 

RPP03990 

RPP04000 

RPP04010 

RPP04020 

RPP04030 

RPP04040 

RPP04050 

RPP04060 

RPP04070 

RPP0408C 

RPP04090 

RPP04100 

RPP04110 

RPP04120 

RPP04130 

RPP04140 

RPP04150 

RPP04160 


cn 

iO 

I 


n o r> 


2 DELTR-0. 

GAHR-TAUP 
GO TO 100 

3 DELTR=S!GNU.tEPSRI*90.*CO^jV 
GAMR=ABSiE?SRl 

GO TO 100 

^ DELTR=S1GN(I..EPSRI*90.«CONV 
GAKR=90.*C0NV-ABSI EPSR > 

100 CONTINUE 
C 

EX^COSIGAMR) 

EV=SIN(GAMRI<‘C£XPI-MINUSJ*DEtTRl 

C 

RETURN 

END 


RPP04170 

RPP041BO 

RPP04190 

RPP04200 

RPP04210 

RPP04220 

RPP04230 

RPP04240 

RPP04250 

RPP04260 

RPP04270 

RPP04280 

RPP04290 

RPP04300 

RPP04310 


SUBROUTINE CLEWEAtEXOf E YOt E XC , EVC, EXXtEYXl 

THIS SUBROUTINE OUTPUTS THE CLEAR WEATHER 1$0LATE0N AND PHASE 

REAL ISOCL 
COMPLEX MINUSJ 

COMFLEX+16 EXOvEYOfEXCsEYCfEXXtEYX 
C0HPLEX*16 VWACtVWAX 
REALMS VWACMfVWAXH 
PEALP8 RTALtAIMAG 
C 

COMMON/BLOC 17M I NUSJ I COvNV 
COMMON/ BL0C2/VU ACM 
C 

VWAC=EXO*EXC + EYO*EYC 
VWAX«EXD*EXX ^ EYO*EYX 
VWACR=REAL(VWACI 



RPP04320 

RPP04330 , 

RPP04340 ^ 

RPP04350 ^ 

RPP04360 , 

RPP04370 

RPP04360 

RPP04390 

RPP04400 

RPP04410 

RPP04420 

RPP04430 

RPP04440 

RPP04450 

RPP04460 

RPP04470 

RPP04480 




c 

c 

c 

*%». 

c 


VMAC1=AIMAGIVUAC) 

VUACM=COABSiVWAC) 

IFdVUACI.EQ.O. }.ANO.( VWACR«EQ.O« ) ) GQ TO I 
VWACPR=ATAN2 I VW ACI tVWACR 1 
VUACPH=VWACPR/CONV 
GO ■^0 2 

1 VWACPH^^O. 

2 CONTINUE 
VUAXR;:=REAt I VHAXI 
VWAXI=AIMAG(VHAX» 

VUAXH-COABSIVWAXt 

IF(VWAX1*EQ.0« .AND. VMAXR.EQ.O.) GQ TO 3 
VWAXPR-ATAN2IVWAX1 >VWAXR» 

VWAXPH=VWAXPR/CDNV 
GO TO 4 

3 Vl^AXPH=0. 

4 CONTINUE 


PHASCW^VWAXPH^VHACPH 


IFIVMAXM.EQ.O. ) GO TO 5 


RPP04490 

RPP04500 

RPP04510 

RPP04520 

RPP04530 

RPP04540 

RPP04550 

RPP04560 

RPP04570 

RPP04560 

RPP04590 

RPP04600 

RPP04610 

RPP04620 

RPP04630 

RPP04640 

RPP04650 

RPP04660 

RPP04670 

RPP04680 

RPP04690 

RPP04700 


I$0CL=20.*DL0G10(VWACM/VWAXMI RPP04710 

MRITE(6,8} ISOCL RPP04720 

8 FORMAT!/, 3X*»THE CLEAR WEATHER ISOLATION REFERENCED TO THE ANTENNARPP04730 
1 PORTS IS* ,19X,F6.2,2X, 'DB* I RPP04740 


GO TO 7 

5 WR[TE(6,6) 

6 FORMAT!/, 3X, 'THE CLEAR WEATHER ISOLATION IS INFINITY* * 

7 CONTINUE 
WR1TE(6,9) PHASCW 


gg 

g £ 
8 1 
fo E 



RPP04750 

RPP04760 

RPP04770 

RPP04780 

RPP04790 

RPP04600 




r, o n o n o n nnoonononri 


9 F&RHAT1/«3X«*THE CLEAR HEATHER PHASE 01 FFEREHCE* >46X« F7.2# 

12X» 'OEGREES*) 

RETURN 

END 

SUBROUTINE ICECGFl FREO.AHOOE ,FXX,FVY»FXy»LGKAN6tRRATE» l.THETAI 

THIS SUBROUTINE RETURNS TO THE CALLING PROGRAM THE EFFECTIVE 
SCATTERING COEFFICIENTS FOR A SLAB OF ICE PARTICLES. THE SLAB 
HAS 50% PROLATE PARTICLES AND 50X OBLATE PARTICLES. THE 
SCATTERING COEFFICIENTS ARE A FUNCTION OF PARTICLE SIZE AND 
ORIENTATION. THE PARTICLES ALL HAVE THE SAME CANTING ANGLE 
THETA. THE COEFFICENTS ARE DERIVED FROM RAYLEIGH SCATTERING 
TECHNIQUES. 

THE PARTICLE ECCENTRICITY IS VERY CLOSE TO ONE. 

DIMENSION ICEflQI.RRATEllSItTHETAf 151 

RIAL LOKANG 

COMPLEX* 16 FXX»FXYfFYY»APROXtAPROYtAOBLXf ADBLYiASPH 


ICE CONCENTRATION INCREASES WITH RAIN RATE. 

IF ONLY CONSIDERING ICE IN THE PATHt RAIN RATE IS MEANINGLESS, 
BUT IT CONTROLS PARTICLE CONCENTRATION. 
RR=1,VOLICE-E-8...RR=50, V0LICE«E-6 

V0LlCE=RRATE(n*{24.4346 08*lE-9l/SQRTIANODE> 

CON*VOL ICE* I FREQ** 21/0.18 
CSTH* COS I THETA C I ) I **2 
SNTH>S1NCTHETA(1 1)**2 
CS L A- COS f L OK ANG I «*2 


RPP04810 

RPP0482Q 

RPP04630 

RPP04840 

RPPO^aSO 

RPP04B60 

RPP04870 

RPP048B0 

RPP04890 

RPP04900 

RPP049I0 

RPP04920 

RPP04930 

RPP04940 

RPP04950 

RPP04960 

RPP04970 

RPP04980 

RPP04990 

RPP05000 

RPP05010 

RPP05020 

RPP05030 

RPP05040 

RPP05050 

RTP05060 

RPP05070 

RPP0S080 

RPP05090 

RPP05100 

RPP051I0 

RPP05120 


I 

•vjf 

fv> 

I 



n o n n r> o o 


5NLA>S1NIL0KANGI*»2 


C 

C 

C 

C 


APAOX= (6 .8 1221 1 ,-0.0266421 
APROY=(3.268519,-0.006179) 

A0BLX=r6.8 12211,-0.026642) 

ACBLY=(2.150047,-0.002674) 

ASPH= 13.954 17,-0. 009044) 

AGBLX=:=ASPH*CSLA *■ A0BLX«SNLA 
A0BLY=ASPH»CSL A ♦ AOBLY^SNLA 

FXX=C0N^{SNTH4J APftOYfrAOBLYl ♦ CSTH«( APRCX^AOBLX)) 
FYY=CON*iCSTH#( APROY+ACBLY) ♦ SNTH4i ( APR QX-fAOBLX ) ) 

FXY=CUN^I APROY*A08LY-APftOX-AOBLX)*SlN(THETACI ) » •COS I THETA f I) J 

RETURN 


RPP05130 

RPP05140 

RPP05150 

RPP05160 

RPP05170 

RPP05180 

RPP05190 

RPP05200 

RPP05210 

RPP05220 

RPP05230 

RPP05240 

riPP05250 

WPP05260 

RPP05270 

RPP05280 


END RPP05290 

SUBROUTINE CAPF( PE , AHOOE ,FXX,F YY.FXY , FREQ, LOKANG,RR ATE ,11, THETA, S I RPPG5300 
IGMA) RPP053I0 


RPP05320 


THIS SUBROUTINE RE TURNS TO THE CALLING PROGRAM THE COMPLEX RPP05330 
SCATTERING COEFFICIENTS FXX,FXY,FYX, FYY. THESE COEFFICIENTS PPP05340 
INCLUDE THE EFFECTS OF A DROP SIZE DISTRIBUTION AND A RPP05350 
CANTING ANGLE DISTRIBUTION. ITHE ABSOLUTE VALUE OF THE CANTING RPP05360 
ANGLE RANGE !S LESS THAN 100 DEGREES) RPP05370 


DIMENSION RR ATE (15), THETA! 15), PE 1 15) ,N(20,28) 
DIMENSION SIGHA(IO) 

INTEGER P,Q 
REAL LOKANG,N 
C0MPLEX*16 FXX,FXV,FYV 
C0MPLFX*16 FV,FH,FSPH,DIFF 


O Q 
« 2 


K V. 


RPP05380 

RPP05390 

RPP05400 

RPP05410 

RPP05420 

RPP05430 

RPP05440 


i 

'-vl 

Co 

I 



n o 


C 

CALL NOIARATEfl I9N^AM00E,THETA»S1G»4A) 

C 

F»X-fO.OfO.OI 

FVY=(0.0f0<^0| 

FXV>(0.0»0^0t 

C 

00 1 P*2, 28 
C 

RAOIUS^O. 125/2.0 ♦ IP-l It'D. 125 
C 

CALL COEF I FREQ. RADIUS. FV«FH,FSPH,0IFF»L0KANG1 
C 


RPP05450 

RPP05460 

RPP05470 

RPP054S0 

P.PP05490 

RPP05500 

RPP05510 

RPP05520 

RPP05530 

RPP05540 

RPP05550 

RPP05560 

RPP05570 


00 1 Q»l,20 

C 

THETAS<^ 1.832594 ♦ 0.174533AQ 
C 

IFfSIGHAf llI.EQ.O.Of THETAS-THETAf 1 1 1 
C 

SNTH*$lNCrHETA$)AA2 
C STH^COS 1 THE TA S > A42 
C 

FXX»FXX ♦ NIQ.PI*! (l.-PEIIII}*FSPHFPEfIlJAfFVASNrH^FH*CSTHn 
FYY*FYY ♦ N(O.PI«f Il-PEI m)AFSPH«^PE(Il )«|FV«CSTH«^FH»SNTHH 
FXV*FXV ♦ Nf0fP)«PE{III«01FFASlN<THETAS}*C0SCTH£TA$l 
C 

1 CONTINUE 
RETURN 
END 

SUBROUTINE NDIRRATE.ll . N .ANODE .THETA .S I6HA I 
THIS SUBROUTINE RETURNS TO THE CALLING PROGRAN 


RPP05580 

RPP05590 

RPP05600 

RPP05610 

RPP05620 

RPP05630 

RPP05640 

RPP05650 

RPP05660 

RPP05670 

RP'^05680 

RPP05690 

RPP05700 

RPP05710 

RPP05720 

RPP05730 

RPP05740 

RPP05750 

RPP05760 



o n o 


C 

C 

c 

c 


c 


c 

c 


c 


THE NUMBER OF DROPS IN THE QTH CANTING ANGLE CLASS AND 
THE PTH DROP SIZE CLASS* 

DIMENSION CANNUNIZOI t AREA! 20) • SIGMAI ?.OI 
DIMENSION RRATEI 15 S jTHETAi 15? ,N120,28I 
INTEGER P,Q 
REAL NtNUMTOT 

NUMTOT=RRATtiIII*5.83333/CAMOOE**3.5l 
1F(S1GMAU1)*EQ«0.0) GO TO 2 
CALL NUMCANINUMTOT aHET A ,CANNUM» U rS IGMA 1 

GO TC A 

2 CONTINUE 
DO 3 J^lt20 
CANNUMI J?=NU«T QT/20.0 

3 CONTINUE 
A CONTINUE 

0ELTAA=0.125 

DO I 0=^1 *20 
DO I P-2.28 

RADIUS=0* 125/2.0 ♦ {P-l)*0.l25 



RPP05770 

RPP05780 

RPP0579G 

RPP05800 

RPP05810 

RPP05820 

RPP05830 

RPP058A0 

RPP05850 

RPP05860 

RPP05870 

RPP05880 

RPP05890 

RPP05900 

RPP05910 

RPP05920 

RPP05930 

RPP059A0 

RPP05950 

RPP05960 

RPP05970 

RPP05980 

RPP05990 

RPP06000 

RPP06010 

RPP06020 

RPP06030 


IF {RADIUS. LE.AMODE ) N( Qt PI=CANNUMl Q)*0ELTAA%RADIUS/AM00E**2 RPP060A0 

IF {RADIUS- GT-AMODE IN *Q . P )= CANNUM* Q )*2.*DELTAA* 1 1 5*R ADI US/AHDDE ) RPP06050 
l/AMODE RPP06060 

IF(RA01US.GT.{2.*AM00E) I N(Q,PI*=0.0 RPP06070 

RPP06060 


I 

cn 

I 



on r* n noor><nn 


! m s 


c 


c 


c 


1 CONTINUE 


RPP06090 

RETURN 


RPP06100 

END 


RPP06110 

SU6R0UTINE NUMCANl NUHT,nT « THETA.CANNUHr 1 1 9 $ IGNA ) 


RPP06120 

RPP06130 

THIS SUBROUTINE RETURNS TO THE CALLING PROGRAH 

THE 

RPP06140 

NUHBER OF DROPS IN ALL THE DISCRETE INTERVALS. 


RPP06150 

THE DISTRIBUTION IS THE NORHAL DISTRIBUTION. 


RPP06I60 

SIGHA CAN EQUAL ZERO OR GREATER THAN 10 DEGREES 


RPP06170 

RPP06180 

OINENSION AREA (201 tCANNUMI 201 


RPP06190 

DIHENSION SIGHAIIOI 


RPP06200 

DIMENSION THETAUSt 


RPP06210 

REAL NUMTOT 


RPP06220 

RPP06230 

SIGMAP^SIGMAI ! n»0.0l7A53 


RPP06240 

RPP06250 

IFISIGHAP.LT.O. 174531 SIGMAP>. 17453 

H 

FC3 0 

RPP06260 


RPP06270 


RPP062B0 

DO 1 I«l,19 


RPP06290 



RPP06300 

K*I-l 

RPP06310 

THETAP--1. 57079 *• K«0. 174533 

& > 

RPP06320 

U-(THETAP-TH£TA«lIU/i«4l42136/SlGMAP 

Is. 

RPP0633G 



RPP06340 

IFfU.LT.O.OI AREA! I »>0.5«ERFC(ABS(U1 1 

*< c.. 

RPP06350 

li^CU.EO.O.OI AREAIII«0.5 


RPP06360 

IFIU.GT.O.QI AREAU)=0.5*ERF(U) * 0.5 


RPP96370 

RPP06360 

1 CONTINUE 


RPP06390 

CANNUHI 1 1 >:NUHTrjT«AREA( 1 1 


RPP06400 


( 

•si 

Oi 

I 


o o r> o r> n n 


C 

C 

c 

c 

c 

c 

c 


c 

c 


2 


CANNUM f 201 =NOMTOTO ( 1 . 0-ARE Af 19 ) ) 
IF(CANNUM(20I.LT*UE-10J CANNUMt 201=0-0 
00 2 1 = 1,18 
K=l + l 

CANNUMLX) = NUMTOT*( AREACKJ-AREAill ) 
lF(CANNUHIK)-LT.l-E-iOi CANNUH(K)=0.0 
CONTIMUF 


RPP06410 
RPP06420 
RPP06A30 
RPP064A0 
RPP06A50 
RPP 06460 
RPP06470 
RPP06480 
RPP06490 
RPP06500 
RPP 06510 
RPP06520 
RPP06530 


RETURN 

END 

SUBROUTINE COEE ( FREO , AMODE ,F V ,EH,F SPH,DIFF ,LOKANGI 

THIS SUBROUTINE RETURNS TO THE CALLING PROGRAM THE SCATTERING 
COEFFICIENTS FOR SPHERICAL AND OBLATE RAIN DROPS. 

THE COEFFICIENTS ARE A FUNCTION OF FREQUENCY AND DROP SIZE 

AND ELEVATION ANGLE. THE COEFFICIENTS USED ARE THOSE OF UZUNOGLU 

EVANS AND HOLT. 

COMPLEX CMPLX 

DOUBLE PRECISION U 1 ,U2 ,U3,U4,U5 
REAL L0KAN6 

C0HPLEX^I6 FV,FH,FSPH,OIFF . 

IFUMODE.lt, 0.25J AM00E=0.25 


RPP06540 

RPP06550 

RPP06560 

RPP065T0 

RPP06580 

RPP06590 

RPP06600 

RPP06610 

,RPP06620 

RPP06630 

RPP06540 

RPP06650 

RPP06660 

RPP06670 

RPP06680 

RPP06690 

RPP06700 

RPP06710 


IF ( IINTI FREQI.EQ.U ) .AND.t AMOOE.GT.3.51 I AH00E=3.5 


RPP06720 


''j 

■vj 



nno nnnnn 




! 


Iff t INTIFREQ),EQ,14).AND,UH00E.GT*3.5) ) AMQ0E*3,5 
IF( ( INT(FREQ)^EQ.20)*AN0,C AH00E.GT.3.0J ) AMODE-3,0 
IFU INT* FREQ).£Q,30I.AND,|AM0DE,GT.3.0) ) AM00E==3.0 
C 

Ul^AMOOE 
U2=AMOOE**2 
U3-AM00E**3 
U4=Af!00E**4 
U5=AMCD£**5 


11,0 GHZ COEFFICIENTS 

IF(IMTlFREQ).NE. 11) GO TO 1 

SPHERICAL DROP COEFFICIENTS 

IFIAHDDE^GTa.OO) GO TO 10 

E0R=-0,0020548 ♦0.01638947*01-0, 041 7568^5J2+0,08832213*U3 
E0I=-0. 0025 154+0, 0192855 3*0 1-0. 04568 16*02+0. 03655 147*03 
GO TO 11 

10 CONTINUE 

E0R=-1. 281 55706+2* 83287 7 18*01-2*07399678*02+0. 60190887*03 
1-0.00984194*04-0.01096165*05 
EOI=2,6027a025-7.52434662*Ul+8.14691632*U2-4. 12133971*03 
1+0.99089467*04-0.08731788*05 

GO TO 12 

11 CONTINUE 



RPP06730 

RPP06740 

RPP06750 

RPP06760 

RPP06770 

RPP06760 

RPP06790 

RPP06B00 

RPP06610 

RPP06820 

RPP06630 

RPP06840 

RPP06350 

RPP06B60 

RPPO607O 

RPP06680 

RPP06890 

RPP06900 ' 

RPP06910 ^ 

RPP06920 

RPP06930 

RPP06940 

RPP06950 

RPP06960 

RPP 06970 

RPP06980 

RPP06990 

RPP07000 

RPP07010 

RPP07020 

RPP07030 

RPP07040 


non 


OBLATE DROP COEFFICIENTS 

EV90R=-0-00 l322+0-01036867*Ul-0. 025372*02+0. 07072533*U3 
EV90I*-0, 0024306+0. 01861967*01-0. 0440732*02+0. 03507413*03 
EH90R--0.0023684+0. 01878307*01-0. 0473704*02+0.09255573*03 
EH90I*-0.0030366+0.02320527*U1-0. 0546296*02+0. 04299093*03 
0IFFR=^0. 00 10464+0. 0084144*01-0. 02 19984*02 +0.0 218304*03 
D IFF 1*-0.000606+0. 0045856*01-0. 0106064*02+0.0079168*03 
C 

GO TO 1000 


RPP07050 

RPP07060 

RPP07070 

RPP07080 

RPP07090 

RPP07100 

RPP07110 

RPP07120 

RPP07130 

RPP07140 

RPP07150 


12 CONTI NOE 
C 

EV90R=-0. 3653892+0. 32540943*01+0.452 80976*02-0. 53907398»'J3 
1+0.20309215*04-0.02432204*05 
EV90I«2.20618555-6. 192 51766*01+6.43873903*02-3.090414*03 
1+0. 7060021 1*04-0.06071 73*05 

EH90R»2. 09022706-7. 07345 1 19*01+9. 036397 26*02-5. 33 156215*03 
1+1.49716423*04-0.15747063*05 

EH90 1=3.28654422-9.682 1 7951*01 + 10. 71 102231*02-5. 5491 184*03 
1+1.36694565*04-0.12512466*05 
D IFF 1=1. 08034369-3.4896 2 08 1*01+4. 2722407*02-2. 45868344*03 
1+0. 66093861*04-0.06440692*05 
C 

IFUHOOE.GT.2.0) GO TO 120 
C 

DIFFR=-0. 32 599997+0. 7531995*01-0.57039996*02+0.15039999*03 
GO 10 1000 
C 

120 CONTINUE 
C 

OIFFR=11.94200 847- 14.59334299*01+5.80600359*02-0.7346671*03 


RPP07160 

RPP07170 

RPP07100 

RPP07190 

RPP07200 

RPP07210 

RPP07220 

RPP07230 

RPP07240 

RPP07250 

RPP07260 

RPP 07270 

RPP0728n 

RPP07?*^ 

RPPO^jO 

RPP07310 

RPP07320 

RPP07330 

RPP07340 

RPP07350 

RPP07360 



r>rir> n o o nnffoorin 


C 

C 


GO TO 1000 
1 CONTINUE 

14 GHZ COEFFICIENTS 
IFI 1NT(FREQI.NE.14) GO TO 2 
SPHERICAL DROP COEEFECIENTS 
IFUMGDE^GT.l.OO) 6C TO 20 

EOR=~0,001376+0.012776*U1-0.04013A*U2+0,128736»U3 
E0i=-0-008796+0*06731467*Ul-0. l58b*U2+0.l24l8133*U3 


GO TO 21 

20 CONTINUE 

EOR^- 12. 13707993*34. 85683676*Ui-37,8408737B*U2+19-48995359*U3 
1-4.69 149894*U4*0.4254B477*U5 

E0I--8. 45501 78*22. ‘79096l7*UI“23.3844400l*U2*ll.09553l6l*U3 
1-2. 41 69453 8*U4*0. 1 9907664*U5 

GO TO 22 

21 CONTINUE 

OBLATE DROP COEFFICIENTS 

EV90R=-0. 000248*0. 0033 5467*UI-O.O13920^U2*O.1OOO2133+U3 
EV90I=-0.008436*0.06455867*UI-0.152256*U2*0.11893333*U3 
EH90R«-0. 000656*0. 007656*Ul-0. 029656«U2*0. I22656*U3 


ROP07370 

RPP07380 

RPP0T390 

RPP07400 

RPP07410 

RPP07420 

RPP07430 

RPP07440 

RPP07450 

RPP07460 

RPP07470 

RPP074B0 

RPP07490 

RPP07500 

RPP075i0 

RPP07520 

RPP07530 I 

RPP07540 -- 

P.PP07550 o 

RPP07560 I 

RPP07570 

RPP07580 

RPP07590 

RPP07600 

RPP07610 

RPP07620 

RPP07630 

RPP07640 

RPP07650 

RPP07660 

RPP07670 

RPP07680 



6H90I^0.010366+0.079154*UI-0«18596*U2*0«143872*U3 
DIfFR=-C^000408+0. 00430133*UI-0. 01 5728*U2*0. 0226346 7*U3 
OIFF I--0^00193*0,0 1459533*01-0. 033 704*U2+0-02493B«.7*U3 

C 

GO TO 1000 
22 CONTINUE 
C 

E V90R =-4.40493663* 12.45481 829*Ul-l2.9018227l*U24^6.32901913»U3 
1-1.4280567 l*U4+0ai87530l*U5 
EV90I=-0. 1923376*0. 3425132*Ul-0. 196242*02*0.0746896 8*U3 
C 

IFtAMf;0E.ST.2.5) GO TO 24 
C 

EH90I=l. 10278458-3. 74696237*U1*4. 4804200 1*U2-2.22964317*U3 
1*0.42010037*04 

D IFFP «-0 , 15268404+0. 39830845*01-0. 3527577*02*0 .11767014*03 
DIFFI— 1.29960041+2.93303624*01-2. 15255009*02*0.5206384*03 
C 

SF2AM00E.GT.2.00} GO TO 25 
C 

EH90R=- 1.02599992+ 2. 30733317*01-1.60799988*02*0.4266664*03 
C 

GO TO 1000 

24 CONTINUE "" 

C 

EH90I=2. 51000033-2.03000022*01+0.62000004*02 

0IFFR=-15. 85152202*18. 27189723*01-6. 68190666*02*0.79411153*03 
OIFF 1=22.3 2321 586-44.48681458*01*30. 17926125*02-8. 46925555*03 
1+0.85300427*04 
C 

GO JO 1000 

25 CONTINUE 


RPP07690 

RPP07700 

RPP07710 

RPP07720 

RPP07730 

RPP07740 

RPP07750 

RPP07760 

RPP07770 

RPP07780 

RPP07790 

RP'»07800 

RPP07810 

RPP07820 

RPP07830 

RPP07840 

RPP07850 

RPP07860 

RPP07870 

RPP07680 

RPP07890 

RPP07900 

RPP07910 

RPP07920 

RPP07930 

RPP07940 

RPP07950 

RPP07960 

RPP07970 

RPP07980 

RPP07990 

RPP08000 


) 


cx> 


I 


r»r>r) n n onr»onnrior> 


r 

C 


EH90R=-19. 65+20. 6567+U 1-6. 66*U2+0.6933*U3 

GO TO 1000 
2 CONTINUE 


\%% 


'll 


20 GHZ COEFFICIENTS 
IFnNTIFREQ).NE.20I GO TO 3 
SPHERICAL DROP COEFFICIENTS 
IFIAMODE.GT.l.OI GO TO 30 

£0R=0. 020296-0. 145276*U1 +0. 297656*U2+0. 008224*03 
eci=-0.0l5488+0. 12709133*U1-0. 334568*U2+0. 30583467*03 




GO TO 31 
30 CONTINUE 

E0R=3. 35567152-8.6 1818659+UI+7. 81 57989* U2-2. 6821 791 7*U3 
1+0.3108073*U4 

E01=»1, 85636463-3, 70750898*Ul+2c 19 599302*U2-0. 25 148 145*U3 


GO TO 32 
31 CONTINUE 

OBLATE COEFFICIENTS 

E V9 OR =0 . 02 4 548 0 0-0 .17912133*01+0, 3 84448* U2-0 • 0 7447467*U3 
CV90I =-0.0 13492+0, 1 1 190533*U1-0.299032*U2+0.27853867*U3 


RPP08010 

RPP08020 

RPP08030 

RPP08040 

RPP08050 

RPP03060 

RPP08070 

RPP08080 

Rppoeooo 

RPP08100 

RPP08L10 

RPP08120 

RPP08130 

RPP08140 

RPP08150 

RPP08160 

RPP08170 

RPP08180 ^ 

RPP08190 00 

RPP08200 

RPP08210 ' 

RPP08220 

RPP08230 

RPP08240 

RPP08250 

RPP08260 

RPP08270 

RPP08280 

RPP08290 

RPP08300 

RPP08310 

RPP08320 



n o r> o o 


EH^0R=0. 02 282-0. 16596667+Ul+O. 34976* U2-0.03061333*U3 
EH90I--0.0 15178+0. 12623933*01-0.33 84 8*U2+0. 31 533867*U3 
DIFFR=0. 000072+0. 00175467*01-0. 015488*02+0. 03426133*03 
DIFFI=-0. 001686+0.014334*0 1-0.039448*02+0.0368*03 


RPP08330 

RPP08340 

RPP08350 

RPP08360 


C 

GO TO 1000 
32 CONTINUE 
C 

EV90R=1. 56037246-3. 97019757*01+3, 66524972*02-1. 24458239*03 
1+0.14526662*04 

EV90I^4. 77376463-13.63610787*01+14.83397997*02 
1-7.66136891*03+1.96881573*04-0. 19749381*05 


RPP08370 

RPP08380 

RPP08390 

RPP08400 

RPP08410 

PrP08420 

RPP0B430 

RPP08440 


C 


C 


C 


EH90R=4. 63026679-1 1.97374109*01+11.0254653*02 
1-4.0050613*03+0.50244373*04 

EH90I--1. 5390057+6. 3417866*01-9.35328879*02+6. 20987895*03 
1-1.74544867*04+0.17820925*05 
DIFFR=3. 06989428-8. 00354339*01 +7. 36021547*02 
1-2.76047886*03+0.3571771*04 

DIFFI=-6. 3 1274803+19.97783334*01-24. 1872053*02+13.8712166*03 
1-3.71425705*04+0.3757024*05 


GO TO 1000 
3 CONTINUE 

If {INHFREQJ.NE.30I GO TO 2000 
30 GHZ COEFFICIENTS 
SPHERICAL DROP COEFFICIENTS 
IF( AMQDE.GT. 1.01 GO TO 40 




I— 

*<, CO- 


RPP08450 

RPP08460 

RPP08470 

RPP08480 1 

RPP03490 -- 

RPP08500 w 

RPP08510 » 

RPP08520 

RPP08530 

RPP08540 

RPP03550 

RPP08560 

RPP08570 

RPP08580 

RPP08590 

PPP08600 

RPP08610 

RPP08620 

RPP08630 

RPP0B64O 



WU (J 


E0R=0. 028004-0. 221564'''U1+0.530744»U2+0* 022816=^ 113 
E01=-0< 007072+0. 09329866*01-0. 390ll2*U2+0. 55688533*03 

GO TO 41 

40 CONTINUE 

E0R=- 1.95096901+3. 4107597*01-0. 82790307+U2-0. 39982^85*03 
1+0.1291964*04 

EO 1=5.77798866-16.3876166*01+16.10102407*02-6.06482662*03 
1+0.82747918*04 

GO TO 42 

41 CONTINUE 

OBLATE DROP COEFFICIENTS 

EV90R=0.03S42- 0.2838733*01+0.7038*02-0. 15834666*03 
EV90I=0*01 3876-0.064756*01-0.026904*02+0.290784*03 
£H90R=0.031652-0.252652*U1+0.614952*U2-0.049952*U3 
EH90 I *-0.000532+0. 04746533*01-0.299512*02+0.50957867*03 
DIFFR=-0. 003 768+0. 03122133*01-0,088848*02+0.10839467*03 
D IFF I*-0. 014408+0. 11222133*01-0.272608*02+0.21879467*03 
C 

60 TO 1000 

42 CONTINUE 
C 

EV90R=3. 93161268-10. 97208846*01+1 1.43870068*02 
1-4. 78953064*03+0.69447892*04 
EV90I=3. 12729055-9. 08249632*01+9. 16246329*02 
1-3.47227982*03+0.47789632*04 
£H90R=-6.97583584+ 16- 79179963*01- 13. 53607686*02 
1+4. 63811854*03-0.57535653*04 


RPP08650 

RPP08660 

RPP08670 

RPP08680 

RPP08690 

RPP08700 

RPP08710 

RPP08720 

RPP08730 

RPP08740 

RPP08750 

RPP08760 

RPP08770 

RPP08780 

RPP08790 

RPP08800 

RPP08810 

RPP08820 

RPP08830 » 

RPP08840 I 

RPP08850 

RPP08860 

RPP08870 

RPP08880 

RPP08890 

RPP08900 

RPPO091O 

ftPPQ8920 

RPP08930 

RPP08940 

RPP08950 

RPP0B960 


EH90I=-1. 39373246+ !• 72736066*U1-0. 14076091*U2+0,05474057*U3 
D IFFR=-XO- 90744851+27. 76388807+UI-24.97477753+U2 
1+9. 4276491 8*U3-1.26983545*U4 
OIFFI-5. 33984552-1 7, 209 20477*Ul+20,78l254l*U2 
1- 11. 64360435»U3+3. 079429 84*U4-0.30370378»U5 


RPP08970 

RPP08980 

RPP08990 

RPP09000 

RPP09010 


C 

1000 CONTINUE 
C 

AL PH= 1 • 5 76 79 6- LOKA NG 
CSLA-0.00l*C0SIALPH)**2 
SNLA-0-001»SIN*ALPHj«*2 
FVR-CSLA*E0R+SNLA*EV90R 
F V I=-CSL A* eO I- SNL A’S'E V90 1 
FHR~CSLA*E0R+SNLA*EH90R 
FHI=“CSLA»E0I-SNLA*EH901 
0IFFR=-SNLA*D1FFR 
DIFFI=SNLA*DIFFI 
EOR-O.O01*E0R 
E0I=-0.001*E0I 
C 


RPP09020 
RPP09030 
RPP09040 
RPP09050 
RPP09060 
RPP09070 
RPP09080 
RPP0‘i090 
RPP09100 
RPP09110 
RPP09120 
RPP09130 
RPP09140 
RPP09150 
RPP09160 
RPP09170 
RPP09180 
RPP09190 
RPP09200 
RPP09210 
RPP09220 
RPP09230 
RPP09240 
ll»14,20,30 GHZ ALLOWED* RPP09250 

RPP09260 

RPP09270 

RPP09280 


FV=JMPLX(FVRfFVI S 
FH==CMPLXIFHR»FHI I 
FSPH^CMPLXCEORiEOl I 
OIFF=sCMPLXCDlFFR,OIFFI ) 

C 

60 TO 3000 
C 

2000 WRITE (6*2001} 

2001 FORM AT I //, 3 X,* FREQUENCY NOT ALLOWED* ONLY 
II 

C 

STOP 



I 

m-t 

00 

Ol 


r»r» o r» ft Oftnnr>nftr> 


3000 

C 


CONTINUE RPP09290 

RPP09300 

RETURN RPP09310 

END RPP09320 

SUDROUTl NE OUTANTI EXN« E YNt EXC t EYC t EXX ■ EVX r I • ISOL 1 • ATTENl « PHASE 1 ) RPP09330 

RPP09340 

THIS SUBROUTINE TAKES THE X AND Y COMPONENTS OF THE WAVE EXITING RPP09360 
THE RAIN CELL IEXN»EVN} AND USES THE COMPLEX VECTOR METHOD RPP09360 

TO COMPUTE VALUES FOR ATTENUATION# ISOLATION# AND PHASE AS A RPP09370 

RESULT OF THE RAIN MEDIUM AND POLARIZATION MISMATCH EFFECTS OF THERPP09380 
RECEIVE ANTENNA RPP09390 

THIS DATA IS STORED IN PROGRAM MEMORY FOR LATER OUTPUT RPP 09400 


RPP09410 

REAL I SOLI RPP09420 

COMMON/BLOCl/MINUSJ.CONV RPP09430 

C0MM0N/BL0C2/VHACH RPP09440 

RPP09450 


DIMENSION PHASElfl 001 # ATTENl 1 100 I # ISOLl I lOOJ 

COMPLEX MI NOS J 
CGMPLEX416 EXC«EYC#£XX#EYX 
COMPLEX* 16 VUPACfVWPAX 
C0HPLEX416 EXNtEYN 
REAL48 VKPACH#VHPAXM 
R£AL«8 REALfAlMAG 

VhPAC-EXN*EXC + EYN*EYC 
VWPAX»EXN*EXX + tVN*EVX 


VHP ACRsR EAL C VHP AC ) 


RPP09460 

RPP09470 

RPP09480 

RPP09490 

RPP09500 

RPP09510 

RPP09520 

RPP09530 

RPP09540 

RPP09550 

RPP09560 

RPP09570 

RPP09580 

RPP09590 


VHPACI^AIHAGIVHPAC ) 


RPP09600 


I 

■ — I 
00 
c» 

J 



VWPACM-COABSI VWPAC ) 

VMPACR-ATAN2 1 V W PAC 1 1 VW PACK ) 

V WPAC P-VWP ACR/COMV 
C 

VHPAXft^REAL IVUPAX) 

VWPAXl:=AlHAGfVHPAXJ 
VWPAXM:=COABS I VWPAX I 
VMPAXR-ATAK2 f V WP AX J « VU PAXR ) 

V WP AX P- V WP A XR/CONV 
1F(VWPAXH.LE.0.000001I VWPAXP^O.O 
C 

PHASE 1 1 VWPAXP-VUPACP 

C 

IF I PHASEU 1 1 .LT. 0. 0) PHASE U II -PHASEl ( 1 1 ^360* 0 
C 

I SOU m ~20.«0L0G1 0 ( VM PACH/ VWPAXH I 
ATTEN 1 1 1 ) -20««0L06 1 0 1 V WACH/ VHP ACM ) 

C 

RETURN 

END 

SUBROUTINE CURVE(X«Yfl 1 I»A«B*R.K) 

C* 

C* THIS SUBROUTINE PERFORMS A LINEAR REGRESSION ON ARRAYS X»Y 
C* 

C* POWER FITS«LOG FITS* AND EXPONENTIAL FITS ARE DONE BY 
C* TAKING THE LOG OF THE INPUT DATA WHERE APPROPRIATE 

C* 

€* THE CALL STATEMENT FOR THIS SUBROUTINE IS 
C* 

C* CALL CURVE(X»Y,Ill»A»BtRfKI 
C* 


RPP09610 
RPP09620 
RPP09630 
RPP09640 
RPP09650 
RPP09660 
RPP 09670 
RPP09660 
RPP09690 
RPP09700 
RPP09710 
RPP09720 
RPP09730 
RPP09740 
RPP09750 
RPP09760 
RPP09770 
RPP09780 
RPP09790 
RPP09800 
RPP09810 
RPP09B20 
RPP09830 
RPP09640 
RPP09850 
RPP09860 
RPP09870 
RPP09880 
RPP09B90 
RPP 09900 
RPP09910 
RPP09920 


> 



I 



c* 

WHERE 



RPP09930 

c* 

X 

(S A SINGLE DIMENSIONED ARRAY WITH 111 

ENTRIES 

RPP09940 

c* 

Y 

IS A SINGLE DIMENSIONED ARRAY WITH III 

ENTRIES 

RPP09950 

c^ 

A 

AND B 

ARE THE RETURNED COEFFICENTS 


RPP09960 

c* 

R 

IS THE 

RETURNED CORRELATION COEFFICENT 

R**2 

RPP09970 

c^ 

K 

IS THE 

FLAG FOR THE TYPE OF CURVE FIT 


RPP 09980 

c* 


K*1 

LINEAR FIT Y=A+B*X 


RPP09990 

c* 


K=2 

POWER FIT Y=A*X**B 


RPP 10000 

c* 


K=3 

LOG FIT Y-A+B*LNIXJ 


RPPlOOlO 

c* 


K=4 

EXP FIT Y-A*EXPIB«XI 


RPP10020 

c* 





RPP10030 


RPP 10040 


C 



* 

RPP10050 


DIMENSION XI3001 »Y(300) 


RPP 10060 

C 



- 

RPP10070 


IF fill. LE. 11 

GOTO 50 


RPP10080 

c 




RPP10090 


Xl~0 



RPPlOlOO 


X2=^0 



RPPlOllO 


V1=0 



RPP10120 


Y2=0 



RPP10130 


XY=0 



RPP10140 


DO 10 l-l.lll 



RPP10150 


XX-XII) 


'&'S. 

RPP10160 


YY=Ym 



RPPI0170 


IFIK.EQ.U 60 

TO 20 


RPP10180 


1FCK«EQ.4) GO 

TO 30 


RPP10190 


XX^ALOGIXXI 



RPP 10200 


IFfK.EQ.3) GO 

TO 20 


RPP10210 

30 

YY=ALOGIYYI 



RPP10220 

20 

Xl^XU^XX 



RPP10230 


X2-X2+XX»KX 



RPP10240 


Yl=Yl+VY 
V2aY2+YY*YY 
XY=XY+XX*YY 
10 CONTINUE 

XNaFLQAimn 

B= I XN* X Y-X 1*YU/IXN*X2-XX^X1) 

AaVl/XN-B*Xl/XN 

Ra ( XN*XY-X l«Yl I ( I XN»X2-XI*X1 ) ^ ( XN*Y2-Y1* Yl I ) 
IF1K.8Q.1.0R.K.EQ.3) RETURN 
AaEXPCAl 
RETURN 
50 A=1.0 

Bal.O 


RPP 10250 
RPP10260 
RPP10270 
RPP10280 
RPP 10290 
RPP10300 
RPP10310 
RPP10320 
RPP10330 
RPP10340 
RPP 10350 
RPP 10360 
RPP10370 


RaQ« 0 

RETURN 

END 


RPP10360 

RPJ>10390 

RPP10400 




I 



